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DETERMINING CENTRAL BLACK HOLE MASSES IN DISTANT ACTIVE GALAXIES. 
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ABSTRACT 

An empirical relationship, of particular interest for studies of high redshift active galactic nuclei 
(AGNs) and quasars, between the masses of their central black-holes and rest-frame ultraviolet (UV) 
parameters measured in single-epoch AGN spectra is presented. This relationship is calibrated to re- 
cently measured reverberation masses of low-redshift AGNs and quasars. An empirical relationship 
between single-epoch rest-frame optical spectrophotometric measurements and the central masses is 
also presented. The UV relationship allows reasonable estimates of the central masses to be made of 
high-redshift AGNs and quasars for which these masses cannot be directly or easily measured by the 
techniques applicable to the lower luminosity, nearby AGNs. The central mass obtained by this method 
can be estimated to within a factor of ~3 for most objects. This is reasonable given the intrinsic un- 
certainty of a factor less than 2 in the primary methods used to measure the central masses of nearby 
inactive and active galaxies, namely resolved gas and stellar kinematics in the underlying host galaxy 
and reverberation-mapping techniques. The UV relationship holds good potential for being a powerful 
tool to study black-hole demographics at high redshift as well as to statistically study the fundamental 
properties of AGNs. The broad line region size — luminosity relationship is key to the calibrations pre- 
sented here. The fact that its intrinsic scatter is also the main source of uncertainty in the calibrations 
stresses the need for better observational constraints to be placed on this relationship. The empirically 
calibrated relationships presented here will be applied to quasar samples in forthcoming work. 

Subject headings: galaxies: active — galaxies: fundamental parameters — galaxies: high-redshift — 
galaxies: Seyfert — quasars: emission lines — ultraviolet: galaxies 



1. INTRODUCTION AND MOTIVATION 

The mass of the central black hole, Mbh, is a funda- 
mental property of active galactic nuclei (hereafter AGNs) 
and quasars governing the physics of their central engine. 
The mass of the central black hole (or "massive dark ob- 
ject") can be measured in nearby, inactive galaxies us- 
ing gas and stellar kinematics (e.g., Kormendy k Rich- 
stone 1995; Richstone et al. 1998; Ho 1999; Kormendy k 
Gebhardt 2001) and in nearby, active galaxies using re- 
sults from reverberation mapping studies (e.g., Koratkar 
k Gaskell 1991; Peterson k Wandel 1999, 2000; Ho 1999; 
Wandel, Peterson, & Malkan 1999; Kaspi et al. 2000). Re- 
cent developments have shown the central mass in both in- 
active and active galaxies to be related to properties of its 
host galaxy, namely the bulge luminosity I*buigei (inactive 
galaxies: e.g., Kormendy k Richstone 1995; Magorrian 
et al. 1998; Richstone et al. 1998; Ho 1999; active galax- 
ies: e.g., Laor 1998; Ho 1999; Wandel 1999; see also Laor 
2001; McLure k Dunlop 2001; Wandel 2001) and the bulge 
stellar velocity dispersion, a (inactive galaxies: Ferrarese 
k Merritt 2000; Gebhardt et al. 2000a; see also Merritt 
k Ferrarese 2001a; active galaxies: Gebhardt et al. 2000b; 
Ferrarese et al. 2001). This is indicative of related forma- 
tion processes of the central mass and the bulge in the 
host galaxy. The reasons for these purely empirical rela- 
tionships, and that the Mbh — & relationship is tighter 
than the Mbh — ^tmigc relationship, are not clear, but 
several attempts to explain them exist (e.g., Silk k Rees 
1998; Haehnelt k Kauffmann 2000; Adams, Graff, k Rich- 
stone 2001; Burkert k Silk 2001; see also § 5 of Merritt 
k Ferrarese 2001b). Further study of these relationships 



and how they may change with redshift is one of many 
motivations for studying supermassive black holes at high 
redshift. 

Spatially resolved gas and stellar kinematical studies 
are most useful for measuring central masses in inactive 
and weakly active galaxies. This is because in AGNs 
the strong glare from the central non-stellar continuum 
source inhibits measurements of the stellar absorption fea- 
tures, and the narrow-line gas kinematics are perturbed by 
non-gravitational forces. Reverberation (or echo-) map- 
ping techniques (e.g., Blandford & McKee 1982; Peter- 
son 2001a) do not require high spatial resolution but in- 
stead utilize the intrinsic continuum source variability and 
light travel time delays between the location of the ioniz- 
ing continuum source and the line-emitting gas respond- 
ing to the changing continuum. If the line-emitting gas is 
gravitationally bound to the central black hole, then sim- 
ple virial arguments imply that the central mass can be 
measured. Reverberation analysis techniques have in re- 
cent years been improved to provide reasonable estimates 
of the size of the broad line region (BLR) and thereby to 
provide reasonable central mass estimates. The accuracy 
of size determinations depend on the data quality, tempo- 
ral sampling, duration of the monitoring campaigns (Col- 
lier, Peterson, k Home 2001; Polidan k Peterson 2001), 
and analyses techniques (e.g., Peterson et al. 1998b; Welsh 
1999; see also Home 2001). Work is in progress to rean- 
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alyze the available AGN Watch data 1 with better analy- 
sis techniques to improve the central mass estimates for 
the Seyfert 1 galaxies which currently have large uncer- 
tainties (B. M. Peterson 2001, private communication). 
Ferrarese et al. (2001) discuss the cause of these large er- 
rors. The simple underlying idea of echo mapping is that 
the measured time delays and the behavior and proper- 
ties of the varying continuum and line emission indeed 
measure distances and velocity dispersion in the BLR and 
thus measure the central mass with reasonable accuracy. 
This assumes that the black hole gravity dominates radi- 
ation pressure effects on the BLR and evidence in favor of 
this now exists (Gaskell 1988; Peterson & Wandel 2000; 
Gebhardt et al. 2000b; Ferrarese et al. 2001). Arguments 
have been presented that it is to be expected, theoretically, 
that other processes than gravity, such as radiation pres- 
sure, anisotropic emission, and projection effects, should 
introduce large systematic uncertainties which dominate 
and prevent accurate measurements of the central masses 
(Krolik 2001). It is worth keeping in mind that theoretical 
considerations are equally affected by our limited knowl- 
edge on the geometry and the details of the kinematics in 
the central regions of AGNs. This seems founded by the 
recent studies, which find good consistencies in the central 
masses of nearby AGNs measured both with reverberation 
and stellar kinematical techniques (Gebhardt et al. 2000b; 
Ferrarese et al. 2001). 

Central masses can now be determined with much 
smaller uncertainties (less than a factor of 2 depend- 
ing 2 on the method and data quality) than the early 
crude estimates could provide (e.g., Dibai 1980; Wandel 
& Yahil 1985; Wandel & Mushotzky 1986; Padovani & 
Rafanelli 1988; Padovani, Burg, & Edelson 1990). The 
accuracy of reverberation mass determinations may fur- 
ther improve given the continuously improving reverber- 
ation analyses techniques. These recent developments al- 
low us to revisit with modern techniques fundamental is- 
sues such as how AGN supermassive black holes evolve 
and relate to galaxy formation and evolution (e.g., Rich- 
stone et al. 1998; Fabian 1999; Kauffmann & Haehnclt 
2000). While reverberation techniques and stellar and gas 
kinematical studies will yield valuable insight to the local 
black-hole demography, the situation is different for high- 
rcdshift AGNs. In principle, the Mbh — c relationship 
is sufficiently tight to allow estimates of the central AGN 
mass based on the measured velocity dispersion, a, of the 
central host galaxy spectrum. This requires, however, a 
significant host galaxy contribution in the AGN spectrum. 
At the time of writing this can only be efficiently studied 
for low- luminosity, nearby AGNs (e.g., Nelson & Whit- 
tle 1995; Ferrarese et al. 2001). Reverberation mapping 
techniques can in principle be applied to distant AGNs. 
Unfortunately, such studies are also very telescope- and 
time-consuming (e.g., Peterson 2001a, b), especially for the 
more distant and fainter AGNs and for the (intrinsically) 



more luminous quasars, which vary on longer time scales 
(e.g., Kaspi 2001). The method is also highly dependent on 
the objects actually showing intrinsic continuum luminos- 
ity variations when observed. As a result, obtaining accu- 
rate reverberation mapping masses of a large, representa- 
tive sample of the distant AGN and quasar populations is 
close to impossible in a human lifetime. Nevertheless, sig- 
nificant advances can be made with reasonable estimates 
of the central black hole masses in these objects even to 
within a factor of a few. This is particularly useful if the 
mass estimates are based on data that are relatively easy 
to obtain such as a single-epoch spectrum. The best avail- 
able approach is likely the "ladder" type of calibration, 
known from the distance scale determination (e.g., Freed- 
man et al. 2001). That is, one calibrates appropriate pa- 
rameters, more easily obtained for higher redshift objects, 
to the mass measurements at low-z. Such an approach is 
pursued here. A relationship between Mbh, the H/3 line 
width, and the optical continuum luminosity, L\(5100A), 
is already established [see eqns. (1) and (2) later] through 
the confirmation of the theoretically expected relationship 
between the BLR size, i?BLR, and L A (5100A) [Wandel 
et al. 1999; Kaspi et al. 2000]. The time delay, t, and 
hence -Rblr = ct for C iv A 1549 is about half that of H (3 
(Korista et al. 1995), and multiple broad emission lines in 
NGC 5548 exhibit the same virial relationship (Peterson 
& Wandel 1999, 2000), as expected if reverberation tech- 
niques work. Therefore, it is reasonable to assume that 
a similar relationship can be obtained with appropriate 
rest-UV spectral measurements. Such a relationship is of 
particular interest if single-epoch spectra can supply the 
required rest-frame UV measurements. Then it becomes 
straightforward to estimate central masses for large sam- 
ples of high-redshift AGNs. If uncertainties in such an 
approach are, relatively speaking, reasonable this is a pow- 
erful tool for statistical studies and black hole demography 
at high redshift. This is the focus of this paper. 

The C iv A 1549 emission line is chosen here to provide a 
characteristic velocity dispersion appropriate for the cali- 
bration of single-epoch UV spectral measurements for sev- 
eral reasons. First, this line is accessible from the ground 
for objects with redshifts between ~1 and ~5. Second, 
its profile is not commonly affected by strong or numer- 
ous absorption lines typical of the Lya emission line 3 . 
Also, its line width is also much less affected by blending 
effects from other line emission, including Fe n emission, 
than either of the C in] A1909 and Siiv+O iv] A1400 lines. 
Contaminating Fe n and He n A1640 emission mainly affect 
the lower profile wings of Civ (e.g., Figure 2 by Marziani 
et al. 1996; Vestergaard & Wilkes 2001; M. Vestergaard et 
al., in preparation). For these reasons it is also an emis- 
sion line on which future rest-frame UV and optical rever- 
beration mapping analyses may be applied in both nearby 
and more distant AGNs, providing tests of the calibrations 
presented here. 



1 All the data from the AGN Watch monitoring campaigns arc available at 
http : //www . astronomy . ohio-state . edu/~agnwatch 

2 This is judged in part from the uncertainties in the masses due to the propagated measurement uncertainties in Rblr an d L\ (5100A) de- 
termined in this work (see also § 7) and the reverberation mass uncertainties quoted by Peterson & Wandel (2000) and Ferrarese et al. (2001). 
Ferrarese & Merritt (2000) quote errors in the stellar kinematics based on the M — <r relationship to be of order 30%, but the scatter can be 
up to a factor of ~2 (Merritt & Ferrarese 2001a). 

3 The exceptions are broad absorption troughs in the blue profile wing seen in a small fraction of quasars. Narrow associated absorption lines, 
detected in some quasars, are generally easily corrected for. 
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It is important to point out that such calibrations do 
not eliminate the need for further reverberation mapping 
studies. Rather, they emphasize the strong need to obtain 
yet better and more monitoring data. First, this work 
and the BLR size - luminosity relationship are based on 
samples of Seyferts and quasars that are, after all, small 
(<34 objects) and do not span the full range of observed 
AGN properties. Second, the UV and optical calibrations 
presented in this work depend strongly on the size — lu- 
minosity relationship and the intrinsic scatter around this 
relationship is the main source of uncertainty in the cali- 
brations (§ 5 and § 6). The presence of this scatter thus 
accentuates the need for more and better BLR size de- 
terminations. It is now quite clear that significant im- 
provements thereof require and can be obtained by a fully 
dedicated set of space-based telescopes, which can provide 
sufficient temporal sampling and duration of the moni- 
toring campaigns (Collier et al. 2001), and simultaneous 
observations across multiple wavelength bands (Peterson 
2001b; Netzer 2001) to constrain the transfer functions 
and the BLR geometry better. In return, more accurate 
measurements can be obtained of the BLR size, and hence 
of the central masses (e.g., Home et al. 2002). This is just 
one of the important issues that the proposed MIDEX mis- 
sion, Kronos, will address if selected. Given the practical 
difficulties of long term multi- wavelength monitoring with 
existing facilities on the ground and in space (HST and 
X-ray telescopes), including obtaining bona-fide simulta- 
neous observing time on many different telescopes for long 
periods of time, Kronos may quite possibly be our best op- 
portunity to study how the central regions of nearby and 
distant AGNs are structured and "engineered" (Blandford 
2001). Ultimately, that will allow more accurate mass es- 
timates for distant AGNs than are possible at present, as 
discussed above. 

This work was inspired by the recently established BLR 
size — luminosity relationship (Kaspi et al. 2000), which 
is applicable not only for nearby Seyfert 1 galaxies but 
also for the intrinsically brighter quasars. Hence, assuming 
this size — luminosity relationship applies to more distant 
quasars also, their central masses can thus be estimated. 
Here, relationships between Mbh and single-epoch rest- 
frame optical and UV measurements, respectively, are de- 
rived and the uncertainties introduced by this approach 
are evaluated. It will be shown that single-epoch spec- 
trophotometry can be used to estimate the central masses 
to within a factor of 3 with high probability, rendering the 
calibrations very useful. 

In what follows 'optical' and 'UV refer to these same 
bands in the rest-frame of the AGN or quasar. A cosmol- 
ogy with H = 75 km s _1 Mpc _1 , q = 0.5, and A = is 
used throughout. 

The paper is structured as follows. The method adopted 
for the calibration is outlined in § 2. Section 3 describes 
the data considered, § 4 discusses how well representative 
the single-epoch optical data arc of the multi-epoch rever- 
beration data, §§5 and 6 present calibrations of single- 
epoch optical and UV measurements, respectively. In § 7 
the performance of the adopted calibrations is evaluated. 



2. THE METHOD 

The virial theorem is used to estimate the mass of the 
central black hole, Mbh ~ rv 2 /G, where v is the velocity 
dispersion of matter at distance, r, which is gravitationally 
bound to the black hole. The velocity dispersion can be 
expressed as v = f ■ v fwhm where i>fwhm is the FWHM 
of the emission profile of the broad line gas. The factor 
/ is on the order of unity and depends on the geometry 
and the details of the kinematics (e.g., Peterson & Wandel 
1999, 2000; Fromerth & Melia 2000; Krolik 2001; McLure 
& Dunlop 2001). The central mass can be expressed as: 



M = 1.5 x 10 6 
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Based on 17 nearby Seyfert 1 galaxies (Wandel 
et al. 1999) and 17 Palomar-Green (hereafter PG; Schmidt 
& Green 1983; Green, Schmidt, & Liebert 1986) quasars, 
Kaspi et al. (2000) determine an empirical relationship be- 
tween the size of the broad line region, -Rblr, and the 
continuum luminosity, AL,\(5100A), where -Rblr is the 
distance of the emission-line clouds responding to the cen- 
tral continuum variations as determined from reverbera- 
tion studies [eq. (6) of Kaspi et al. 2000] : 
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In appendix A this relationship is re-derived using re- 
gression analysis appropriate for data with internal scat- 
ter. A critical assessment of the object sample is also per- 
formed. The modified luminosity dependence has a larger 
uncertainty, which is probably more realistic given the in- 
trinsic scatter, than that quoted above. However, within 
this uncertainty the approach adopted in appendix A does 
not yield an -Rblr — L relationship significantly different 
from that of Kaspi et al. for the current sample of AGNs; 
this may change once more data are available to better 
constrain the -Rblr — L relationship. Therefore, in what 
follows cqn. (2) is used to estimate the BLR size when 
unknown. 

Equation (1) assumes that the broad line emitting gas is 
gravitationally bound, and that the cloud velocity disper- 
sion is isotropic such that v = \/3|<7j|, i = 1,2,3; |<Tj| = 
wfwhm/2; wfwhm is in units of kms -1 . The details of the 
geometry and kinematics, i.e., the exact value of /, are 
unknown but have been debated (e.g., Peterson & Wandel 
1999; Fromerth & Melia 2000; McLure & Dunlop 2001). 
Most expected values of / are anticipated to result in a 
constant offset in log Mbh (Krolik 2001). Because such 
details are not yet clarified or well constrained (e.g., Wan- 
del 2001) the approach taken in this work is to compute 
the central masses with basic and common assumptions 
to allow comparison with most other work. The offsets 
in mass estimates due to different assumptions of geome- 
try and kinematics should however be kept in mind when 
considering the Mbh values in absolute terms. 

Assuming that equation (2) is valid for all active galax- 
ies, we therefore have an approximative relationship be- 
tween the mass of the central black hole, the AGN con- 
tinuum luminosity at 5 100 A and the line width of the H/3 
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emission-line [eqn. (1)]. Given that the BLR size — lu- 
minosity relationship is determined from line-continuum 
variability of the objects, the specific H/3 line width to be 
used to determine the 'reverberation mass' is that appro- 
priate for the emission-line gas at the distance -Rblr from 
the ionizing continuum source that is varying. That is, 
the FWHM of H (3 in the 'root-mean-square 4 (rms) spec- 
trum' (e.g., Peterson & Wandel 1999, 2000) will be used 
below for the variability data. Kaspi et al. (2000) argue 
that FWHM(H/3, mean), the FWHM of H/3 in the 'mean 
spectrum', may equally well be used. In the interest of 
completeness and for comparison, the masses based on 
the mean optical spectra are included in the analysis be- 
low, even if the rms masses are considered, strictly speak- 
ing, the most representative of the intrinsic, actual central 
masses. 

The method used here to calibrate single-epoch UV 
measurements to reflect the central black-hole masses is 
briefly as follows. First, the single-epoch FWHM(H/3) 
and continuum luminosities are compared with the multi- 
epoch equivalents to test how representative the single- 
epoch measurements are, and whether systematic offsets 
between these two sets of measurements exist (§ 4). If 
that is the case, a calibration or modification of the single- 
epoch measurements may be appropriate to ensure they 
introduce minimum systematic error or do not add un- 
necessary scatter. Then, the single-epoch mass estimates 
based on optical measurements, Mbh(H/3), are computed 
using equations (1) and (2). These estimates are com- 
pared and calibrated to the available reverberation masses, 
-^bh (H/3, rms), based on H/3 measurements for the PG 
quasars (§ 5). The next step is to determine the best 
calibration of an appropriate combination of the single- 
epoch FWHM(Civ) and A L A (1350A) measurements to 
the best available Mbh measurements (§6). This is done 
for the subset of 26 objects with both MBH(H/3,rms) and 
UV measurements available. How well this method yields 
representative central mass estimates is discussed in § 7. 

3. DATA 

The calibration of optical single-epoch spectral measure- 
ments is based on the 19 objects common to the sam- 
ple of Seyfcrt Is and PG quasars with established rever- 
beration masses (Peterson & Wandel 1999, 2000; Wandel 
et al. 1999; Kaspi et al. 2000) and the PG quasars studied 
by Boroson & Green (1992; hereafter BG92). The cali- 
bration of the single-epoch UV measurements is based on 
the sample (labeled "UVrev") of 26 reverberation AGNs 
and PG quasars (Wandel et al. 1999; Kaspi et al. 2000) 
with available UV spectroscopy. Mostly for comparison, 
a sample is also considered consisting of UVrev and 30 
additional PG quasars (BG92) for which C iv line widths 
and UV continuum luminosities are also available in the 
literature. These data are described in the following. 



the contaminating iron emission and the narrow line con- 
tribution are subtracted; i.e. the FWHM is that of the 
supposedly intrinsically emitted broad H/3 line emission. 
Some of the line widths were substituted with measure- 
ments obtained in this work, because the subtraction of 
the narrow component renders the profile less representa- 
tive of the line width needed for this study (§ 4.1). 

Monochromatic luminosities, La(5100A), are deter- 
mined for most of the PG quasars from the spectropho- 
tometry by Ncugebauer et al. (1987). The flux densities 
at 5100A are approximated by a linear interpolation. For 
most objects this interpolation is done over short or neigh- 
boring frequency ranges in relatively tight spectral en- 
ergy distributions, and so should be reliable. For the few 
sources not included in this database La(5100A) is com- 
puted based on the 4400A flux densities listed by Keller- 
man et al. (1989). These are based on B-band photome- 
try of Schmidt & Green (1983). A power-law continuum, 
F„ ex y Q , with an average optical slope, a — a opt — —0.5 
is assumed and extrapolated to 5100 A. The optical data 
are listed in Table 1. 

Errors on the parameters are propagated using the 
following assumed uncertainties: (1) spectral slopes, 
cr(ai)=0.2, reflecting the different slopes commonly 
adopted in the literature for quasars (e.g., Fran- 
cis et al. 1991: a=-0.3; Veron-Cetty & Veron 1993: 
a=— 0.7), (2) redshift, cr(z)=0.0025, a compromise be- 
tween the typical uncertainty for redshift determinations 
based on optical and low-ionization lines, and those relying 
on UV high-ionization lines only where velocity shifts may 
exist between lines, especially for higher redshift quasars 
[cr(z) has a negligible effect], (3) Neugebauer et al. (1987) 
quote flux errors for each measurement and for each object, 
and (4) the B-magnitude uncertainty is 0.27 mag (Schmidt 
& Green 1983). BG92 do not quote the uncertainty on 
their FWHM(H/3) measurements. Based on the few error 
quotes available in the literature (e.g. Brothcrton 1996; 
Vestergaard 2000; M. Vestergaard et al. , in preparation) 
a reasonable relative error on the FWHM measurements 
of broad emission lines is ^10% depending on the mea- 
surement method and the quality of the data. Blending 
effects will increase this uncertainty. The uncertainties in 
the BLR size — luminosity relationship [eqn. (2)] are also 
taken into account; see also Appendix A. 

3.2. UV Data 

The single-epoch UV measurements were obtained from 
a handful of studies including Wilkes et al. (1999), Wang, 
Lu, & Zhou (1998), Laor et al. (1994, 1995), Koratkar & 
Gaskell (1991), and Wills et al. (1995). Line widths and 
luminosities were generally taken from the same study. 
For a few objects with data only presented by Marziani 
et al. (1996) the line widths were remeasured directly on 
the published C iv profiles. This was done because the au- 
thors decompose the profiles and only list measurements 
for the individual broad and narrow components. As ar- 
gued in § 4.1 the line widths for this study should be 
measured on profiles including the narrow component as 

4 Based on the spectral database obtained during a monitoring campaign a 'mean spectrum' can be generated, which is the mean flux at a 
given wavelength bin of all the individual spectra. Similarly, the root-mean-square deviation at each wavelength bin from this mean can be 
obtained. This rms spectrum shows how the object spectrum varies in response to continuum emission variations. See Peterson et al. (1998a) 
for details and examples of 'rms' and 'mean' profiles. 



3.1. Optical Data 

BG92 list H/3 FWHM measurements of single-epoch 
spectra of many quasars in the PG sample from which 



Black Hole Masses in Distant Active Galaxies. 



5 



it originates in the BLR like the broad component; this 
is contrary to the assumption of Marziani et al. (1996). 
Most of these earlier studies present UV continuum lumi- 
nosities at 1350A or quite close to it. Wilkes ct al. (1999) 
quote continuum luminosities at 2500 A, extrapolated from 
available V and B band photometry. Most of the objects 
studied by Wilkes et al. were also studied by others who 
do measure the continuum luminosities directly. These 
latter measurements are less subject to extrapolation un- 
certainties and were thus preferred. All the luminosities 
adopted from the literature were extrapolated to be a mea- 
sure of the 1350 A luminosity, if necessary, assuming the 
slope and cosmology quoted above; this extrapolation is in 
all cases over a relatively short wavelength range (<100A). 
For the measurements with quoted uncertainties the er- 
rors were propagated, as outlined above. For line widths 
with no measurement errors a conservative uncertainty of 
10% was imposed. Typical luminosity uncertainties of 30% 
were adopted for luminosities without quoted errors. The 
discussion in § 4.2 show this is reasonable. 

For some of the objects two or more studies quote line 
widths and luminosities. It was then assumed that the dif- 
fering measurements is a consequence of intrinsic scatter 
due both to variability and to measurement uncertainty. 
The mean of the measurements was therefore adopted and 
the difference in line width was adopted as the uncertainty 
in the FWHM measurement, unless it was smaller than the 
adopted conservative 10% error. The luminosities were 
generally rather similar. The adopted UV measurements 
and their source references are listed in Table 2. 

4. CALIBRATION OF SINGLE-EPOCH OPTICAL 
MEASUREMENTS 

Before the single-epoch spectral measurements (line 
width and monochromatic luminosity) are calibrated to 
the (multi-epoch) reverberation measurements it is im- 
portant to compare these two types of measurements to 
examine whether single-epoch measurements are reason- 
ably representative of the reverberation measurements, to 
what degree and whether certain spectral corrections are 
necessary. This includes an attempt to understand how 
and in what way they may differ. This direct compar- 
ison is performed on the objects common to the Kaspi 
et al. (2000) and BG92 studies, namely the 17 PG quasars 
studied by Kaspi et al. and two Seyfert galaxies, MrkllO 
and Mrk335, (Wandel et al. 1999), all of which were stud- 
ied with reverberation mapping techniques. 

4.1. Line Width Measurements 

The appropriate BLR velocity dispersion to use in the 
mass estimate is that measured from the 'rms' profile 
(§ 2). Because the mass depends strongly on the line 
width [cqn. 1] it is important to identify the single-epoch 
line profile that best represents the rms profile. That is, 
is it important to subtract contaminating Fe n emission, 
and/or perhaps the narrow line contribution before mea- 
suring the FWHM? Another goal of this section is to quan- 
tify how well a single-epoch line width measurement esti- 
mates the rms line width. Rms spectra show the variable 
part of the continuum and line emission and its strength 
as a function of wavelength. The H/3 line width in the 
rms spectrum, hereafter FWHM(H/3, rms), is the width 



of the variable part of the line. Rms spectra can have dif- 
ferent line emission contribution from single-epoch spec- 
tra (examples of mean and rms profiles are presented by 
Peterson et al. 1998a). The H/3 line width measured in a 
single-epoch spectrum, hereafter FWHM(H /3), is expected 
to be more similar to the mean of many individual (single- 
epoch) spectra than FWHM (H/3, rms), which depends on 
which part of the broad line gas that varies. Therefore, in 
what follows, after a general discussion of the single-epoch 
and multi-epoch line widths, the FWHM(H /3) values are 
first discussed in relation to FWHM(H/3, mean) and then 
to FWHM(H/3, rms). 

4.1.1. The Best Single-Epoch Profile to Measure 

Figures la and lb show the FWHM(H /3, mean) and 
FWHM(H/3, rms), respectively, and their uncertainties 
measured by Kaspi et al. and Wandel et al. in their 
multi-epoch spectra are plotted for the 18 objects of 
those mentioned above with H/3 measurements versus 
the single-epoch FWHM(H/3) measurements by BG92. 
Note that there is a typographical error in Table 2 of 
BG92, FWHM(H/3) = 5320 kms" 1 for PG1307+085 (Laor 
2000). The dotted line represents a one-to-one relation- 
ship. With a few exceptions (labeled in Fig. 1) the 
single-epoch FWHM(H/3) shows a general consistency 
with FWHM(H/3, mean) within 15% - 20% variation in 
the single-epoch width and also show consistency with 
FWHM(H (3, rms) to within -20% - 25% variation. This 
is consistent with the ±15% line width variation (B. 
M. Peterson 2001, private communication) observed for 
NGC5548 during 1988 December 14 - 1996 October 16 
(e.g., Peterson et al. 1999), where the H/3 line flux varied 
up to ±66%. It indicates that most of the single-epoch 
FWHM(H /3) measurements and their deviation from the 
multi-epoch measurements can perhaps be accounted for 
as being due to intrinsic continuum and line variations. A 
detailed discussion of these deviations follows below. 

It is important to note that there are technical differ- 
ences between the BG92 and the multi-epoch measure- 
ments. BG92 fit and subtract the optical Fen emission 
around H/3 and [Oin] AA4959,5007 and also fit and sub- 
tract the narrow line H /3 contribution. Kaspi et al. and 
Wandel et al. do not perform any of these corrections to 
their data. Subtraction of the narrow line contribution 
will increase the FWHM value. This increase may be quite 
significant depending on the strength of this contribution 
(see below and e.g., Jackson & Browne 1991). By elimi- 
nating the Fen emission, which blends into the red wing 
of H/3 and both of the [Oin] doublet lines, a smaller H/3 
width is obtained. This holds as long as the contaminating 
Fe II emission does not artificially increase the underlying 
continuum level significantly, creating a so-called "pseudo- 
continuum" (e.g., Wills, Netzer, & Wills 1985; Vestergaard 
& Wilkes 2001). Subtracting the Fell emission will then 
tend to increase the FWHM as the line peak height is in- 
creased (the local continuum level or "zero-point" of the 
line is lowered). Thus, the effects of subtracting the Fen 
emission is not straightforward to predict a priori. 

For this type of study the FWHM(H/3) measurement 
that best represents FWHM (H/3, rms) is one measured on 
a profile for which only a correction for strong Fell con- 
tamination is performed, and especially one in which the 
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narrow line contribution is not subtracted. This is justified 
in the following. Note that if the Fen emission is essen- 
tially constant in strength while the object was monitored, 
the Fe n emission contribution is automatically eliminated 
in the rms spectra. AGN optical Fen emission does not 
seem to vary much or very fast if at all (e.g., Wamsteker 
et al. 1990; Kollatschny & Welsh 2001) so it is reasonable 
to assume that the FWHM(H/3, rms) is not significantly 
affected by Fell emission, if present. For this reason it is 
desirable to use single-epoch FWHM(H/3) measurements 
which are not affected by the presence of Fe n emission. 

For these reasons and because some of the BG92 
FWHM(H/3) measurements deviate strongly from the 
FWHM(H /3,rms) of Kaspi et al. and Wandcl ct al. (see 
below and Fig. 1), the line width was remeasured in this 
work in the original spectra (i.e., those not corrected for 
Fe II or narrow component emission; these data were kindly 
provided by T. Boroson) of the PG quasars studied by 
Kaspi et al. With exception of a handful of objects with 
strong Fen emission and/or a strong narrow line emis- 
sion component, these non-corrected FWHM values are 
consistent with the BG92 values to within the (assumed) 
10% measurement uncertainties; the inconsistent measure- 
ments were corrected as described below. The objects 
showing "uncorrected" FWHM(H (3) values larger than the 
BG92 measurements are all Fen-strong, indicating that 
the Fe II blending effects in the red wing of H j3 dominate 
the effects of the alleged Fe n-pseudo-continuum for these 
objects. For these objects the Fen corrected BG92 mea- 
surements were chosen to be the most representative of 
FWHM(H (3, rms) since this width is not significantly af- 
fected by Fe n emission. However, the uncertainties of the 
BG92 FWHM(H (3) measurements for these objects were 
corrected to reflect the effect of the Fe n emission. The 
error was set to the difference between the measurements 
before and after the Fe n emission (and narrow component) 
subtraction if this difference was larger than the conserva- 
tively assigned 10% error. The most strongly Fen con- 
taminated quasar, PG1700+518 show an FWHM of ~650 
kms -1 larger than the BG92 value (Fig. 1). 

PG1704+608 clearly illustrates why the narrow line core 
should not be subtracted. PG1704+608 has a very strong 
narrow component and exhibits a highly significant offset 
in FWHM(H/3) from FWHM(H/3, mean). Its rms spec- 
trum shows that the strongest H (3 variation furthermore 
occur in the narrow line core [the FWHM(H (3, rms) is 
quite small; Figure lb], explaining why the BG92 FWHM 
is far from representative of the typical FWHM(H/3, 
mean) and FWHM(H/3, rms) measured for this object. 
This means that the low- velocity gas emitting H (3 is in- 
deed well within the BLR and that the narrow line core 
is bona-ftde BLR line emission. Typical variability time 
scales of the narrow emission component support this 
(Stirpc 1990). For PG1704+608 the BG92 FWHM(H f3) 
was therefore replaced, in the analysis described below, by 
the FWHM(H/3) measured here in the uncorrected spec- 
trum. 

4.1.2. Line Width Deviations 

The objects labeled in Figure la (except PG1617+175) 
have BG92 FWHM(H/3) measurements larger than 
FWHM(H/3, mean) in excess of ^20% variation and the 



typically expected measurement uncertainties. Are these 
deviating measurements easily understood? PG0052+250 
and PG1307+085 are quite variable in the H/3 line flux 
as indicated by the light curves of Kaspi et al. (2000) and 
so their deviating line widths are not unexpected. The 
"uncorrected" FWHM(H (3) of MrkllO and PG1426+015 
deviate from the BG92 measurements by more than 10%. 
Both objects have strong (spiky) narrow H f3 components 
(see e.g., Fig. 1 by BG92) that were subtracted by BG92. 
This explains the larger BG92 widths. PG1613+658 be- 
haves strangely. The mean H (3 profile is broad but it varies 
most strongly in the narrow line core (see e.g., Figure lb 
and Table 6 and Figure 5 of Kaspi et al. 2000). However, 
for this object the narrow component of neither the single- 
epoch profile nor the mean profile is a fair representation 
of the rms profile (Figure lb). Therefore, this object will 
commonly be excluded from the analysis. To reiterate, for 
Mrk 110, PG1426+015, and PG1704+608 the single-epoch 
FWHMs based on the original, uncorrected data are used 
in the analysis instead of the BG92 values. For the re- 
maining objects the single-epoch FWHM(H/3) is consis- 
tent with FWHM(H (3, mean) to within reasonable effects 
of line variability. 

Figure lb shows the single-epoch FWHM(H (3) of BG92 
with the corrections discussed above plotted versus the 
FWHM(H/3, rms) [Wandel et al. 1999; Kaspi et al. 2000]. 
Most of the objects with single-epoch FWHM(H (3) deviat- 
ing from FWHM(H/3, rms) by more than +15% were well 
within ±15% from the FWHM(H/3, mean) in Figure la. 
This may indicate that the variation in FWHM is typically 
a little larger than 15%, as is seen for NGC 5548, in these 
objects and mostly occur in the narrow line core. The 
main technical difference between the BG92 FWHM(H f3) 
and the FWHM(H/3, rms) is that the BG92 FWHM only 
reflects the width of the broad line component as both line 
widths are expected to be little affected by the presence 
of Fe II emission. For objects with strong variation in the 
low-velocity gas the two FWHM measurements will there- 
fore differ significantly as observed and discussed above. 
Of the two objects (PG1617+175 and PG1613+658) with 
single-epoch FWHM(H/3) (and error) deviating by >15% 
from FWHM(H/3, rms), only PG1613+658 has 'such a 
significantly deviant BG92 FWHM(H/3) that the offset 
cannot be ascribed to larger measurement uncertainties 
or stronger source variation than the ±15% shown by 
NGC 5548. 

It is worth noting for completeness that the differences 
in spectral resolution between the BG92 and the Kaspi 
et al. spectra do not contribute to the line width differ- 
ences. The original BG92 spectra were degraded to the 
resolution of ^10 A of the Kaspi et al. spectra for this 
analysis. No significant FWHM differences were found. 

4.1.3. Regression Analysis 

Regression analyses were performed to statistically test 
how well the single-epoch FWHM(H/3) represents the 
multi-epoch line widths. It is clear from Figure la that the 
objects with line widths below 4000 kms -1 yield very con- 
sistent measurements in spite of their variability, while the 
broader lined objects have larger uncertainties from both 
measurements [recall, that a 10% measurement error was 
assumed for the single-epoch FWHM (H/3)] and especially 
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from intrinsic source variability; these objects have rather 
similar luminosities to the other objects in the sample, 
and so these differences are not likely luminosity-related. 
However, PG1613+658, PG1426+015, and PG1617+175 
have the broadest mean H line and also strong, broad 
Fen emission (see BG92). This emphasizes the need for 
linear regression methods, which are more appropriate for 
the nature of these data. The regression analyses were 
therefore performed using the bivariate correlated errors 
and intrinsic scatter (hereafter BCES) algorithm (Akritas 
& Bershady 1996). When intercomparing the multi-epoch 
measurements with those of single-epoch spectra, intrinsic 
scatter and measurement errors are expected. The BCES 
algorithm is the most appropriate to use as it does not, as 
do many other linear regression methods, assume a per- 
fect relationship between the variables if the measurement 
errors could be made insignificantly small. Mcrritt & Fer- 
rarese (2001a) compare the BCES algorithm to other com- 
monly used linear regression methods. 

The uncertainties used in the regressions are the sym- 
metric errors on the logarithms based on the positive linear 
errors. No significant differences were found when using 
the negative linear errors instead. These symmetric er- 
rors [= (log e) (ax I 'X)] were determined by propagating 
the linear errors to the logarithmic relationship. 

The best fit to the FWHM(H 0, mean) and FWHM(H 0) 
distribution is shown in Table 3. It is based on sample B, 
that is, the PG quasars (Kaspi et al. 2000) and MrkllO 
and Mrk335 (Wandel et al. 1999), excluding PG1351+640 
as it has no H/3 measurements available. The BCES re- 
gression is quite robust; bootstrapping simulations repro- 
duce the theoretically expected results well (see e.g., Akri- 
tas & Bershady 1996). Also, there is little difference be- 
tween the BCES(Y|X) [i.e., Y = f(X)] and BCES(X|Y) 
[i.e., X = g(Y)] regressions. The bisector bisects these two 
regressions and so the said small difference is reflected in 
the similar results obtained for the BCES(Y|X) and bisec- 
tor regressions (both are listed in Table 3). These BCES 
regressions for all the data points in Figure la (sample 
B) are consistent with a unity relationship within 1.8c 
When the FWHM(H 0, rms) measurements are considered 
(Fig. lb) PG1613+658 is very much an outlier and it is 
excluded from the regressions to those data as it behaves 
in an obviously strange fashion (§ 4.1.2). For comparison 
with the results of that analysis, the regressions were re- 
peated on the sample in Figure la excluding this object 
(sample C): the data now display a tighter relationship and 
the BCES regression is consistent with a unity relationship 
to within lc, with no regard to which is the independent 
variable (Table 3). Only the BCES bisector for sample C 
is shown in Figure la for visibility. The BCES bisector 
is also the most representative of the intrinsic relationship 
between the two line widths, because they are different 
measures of the same property. 

In Figure lb the bisector regression (solid line) for 
sample C shows consistency with a unity relationship to 
within lc FWHM(H/3, rms) = (0.98±0.08)x single- 
epoch FWHM(H 0) - (265±227). Whether this relation- 
ship or a pure 1:1 relationship is used, the propagated 
effect on logMeH is only of order ^0.04 dex for a 4000 
kms" 1 line width (and ranging between 0.07 and 0.03 dex 
for widths of 2000 kms -1 and 6000 kms -1 , respectively) 



which is well within typical FWHM measurement uncer- 
tainties. So, assuming a 1:1 relationship is fair, as long 
as the single-epoch line width is measured as described in 
§ 4.1.2. 

4.1.4. Is the High-Velocity BLR Gas Optically Thin? 

It is interesting to note that this comparison of 
FWHM(H/3, rms) and single-epoch FWHM(H/3) (with 
or without the Fen emission) shows that for the PG 
quasars the single-epoch FWHM(H 0) is predominantly 
larger than the FWHM(H 0, rms). If all the line flux, both 
at high and low velocity, in the profile is varying with simi- 
lar amplitudes one would naively expect the rms profile to 
be consistently broader than any given single-epoch pro- 
file as the rms profile represents the responding amplitudes 
and the velocities of the responding gas. At the very least, 
the single-epoch width should scatter around FWHM(H 0, 
rms). Yet the opposite is seen here. One may speculate 
whether this indicates that the broad line wings, which are 
often not represented in the rms line profiles, are mostly 
due to optically thin material, which does not respond to 
continuum variations as does optically thick gas? This is 
probably seen for He II in NGC 4051 (Peterson et al. 2000). 
PG1704+608 and, especially, PG1613+658 were pointed 
out above as extreme cases showing broad average pro- 
files, yet exhibit variability in the narrow line core only. 
Other objects showing similar behavior in their H line 
profiles are Mrk590 (Fcrland, Korista, & Peterson 1990) 
and Mrk335 (Kassebaum et al. 1997). Similar behavior 
has been seen for both the Ly a and C iv line profiles in the 
higher luminosity quasars (e.g., O'Brien, Zheng, & Wilson 
1989; Perez, Penston, & Moles 1989; Gondhalekar 1990). 
Therefore, the presence of high- velocity, optically thin line 
emission is likely rather common in AGNs and quasars. 
Shields, Ferland, & Peterson (1995) discuss this issue and 
its possible importance. They further argue that the pres- 
ence of such optically thin emission gas can explain some 
of the variability properties of Seyfert Is. 

To summarize § 4.1, statistically, the single-epoch 
H.0 line width displays a unity relationship with both 
FWHM(H/3, mean) and FWHM(H 0, rms) to within 
lc However, this is valid as long as the single-epoch 
FWHM(H/3) is measured on H/3 profiles, which includes 
the narrow component but for which especially strong (and 
blending) Fell emission is subtracted. 

4.2. Luminosity Measurements 

The optical continuum luminosities for the single-epoch 
spectra are gathered mostly from Neugebauer et al. (1987) 
supplemented with data from Schmidt & Green (1983) for 
Mrk 110 and Mrk 335, as described in § 3. In Figure 2 these 
ALa(5100A) measurements are compared to those deter- 
mined by Wandel et al. (1999) and Kaspi et al. (2000) for 
the same sample of quasars discussed above in § 4.1 but 
including also PG1351+640 (i.e, sample A; Table 3); the 
luminosities are plotted with the same cosmology (Ho = 
75 km s _1 Mpc _1 , qo = 0.5, and A = 0). The uncertainties 
listed by Kaspi et al., based on the rms in the continuum 
light curves, are also plotted. The single-epoch luminos- 
ity uncertainty reflects the propagated errors based on the 
measurement uncertainties as described in § 3.1. It is clear 
from Figure 2 that for most of the objects the luminosity 
measured by Neugebauer et al. is stronger. 
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Regression analyses show (Table 3) a best fit slope con- 
sistent with 1.0 and an intercept of zero to within the 
uncertainties. The various BCES results are rather robust 
and show little difference. Figure 2 indicates that most of 
the data points are scattered around a slope of one with 
a systematic offset of ~0.1 dcx, (i.e., <~30%). This offset 
is equivalent to an offset in logM of 0.07 dex. As will be 
clear later, this is well within the uncertainties of the (cal- 
ibrated) mass estimates and will thus not significantly af- 
fect the calibration of optical single-epoch mass estimates. 
Therefore the central mass estimate calibration is eval- 
uated assuming no luminosity offset, i.e., no correction 
of either luminosity measurements is applied, but keeping 
the offset in mind. The ~0.1 dex scatter in the luminosity 
measurements around the regression line is adopted as a 
representative uncertainty in using single-epoch luminos- 
ity measurements to measure the AGN mean luminosity 
(thus justifying the choice in § 3). A higher uncertainty 
may apply in practise, especially if (spectro-)photometric 
data are not used to determine the continuum luminosity. 

For completeness it is noted that the luminosity off- 
set is not caused by either (1) interpolation errors in 
the Neugebauer et al. data [tight spectral energy distri- 
butions; interpolation across neighboring pixels], (2) in- 
trinsic source variability [possibly some contribution, but 
symmetric scatter around a one-to-one relationship would 
be expected], (3) aperture differences between the Kaspi 
et al. and Neugebauer et al. studies, or (4) different cor- 
rections for reddening or Galactic extinction. As for (3), 
both studies have apertures including equal amounts of 
host galaxy contributions, which afterall is only strong in a 
few objects (S. Kaspi 2001, private communication). How- 
ever, the use of slightly different absolute flux calibration 
scales could be the origin of the systematic luminosity off- 
set. 

5. CALIBRATION OF SINGLE-EPOCH MASS ESTIMATES: 
OPTICAL MEASUREMENTS 

Having established that the single-epoch line widths and 
luminosities are statistically representative of the multi- 
epoch equivalent measurements, the "single-epoch central 
mass estimates" can be computed and further analyzed. 
The first step in the calibration of UV single-epoch mea- 
surements to estimate the central AGN/quasar mass is to 
compare and calibrate, if necessary, the single-epoch opti- 
cal measurements to yield representative estimates of the 
central masses. This is done: here using the sample: of 18 
objects with central mass determinations from reverbera- 
tion mapping (Wandel et al. 1999; Kaspi et al. 2000) and 
H measurements. 

Using the line widths and luminosities as discussed 
above the single-epoch mass estimates, Mbh(H/3), were 
computed for the 18 AGNs and quasars using equations (1) 
and (2). These masses are plotted in Figure 3 versus the 
masses determined from reverberation mapping studies. 
The masses quoted by Kaspi et al. (2000) are based on 
the average of the H/3 and Ha line widths. It is more 
appropriate to compare the single-epoch mass estimates 
with reverberation masses determined from the H line 
width only. The reasons are the following: (1) Kaspi 
et al. determined the BLR size, Rblr, using data on both 
H a and H but found the same results using H alone, 



just with larger scatter, (2) The H a A 6563 profile is poten- 
tially contaminated by [N n] A A 6548, 6583 which will affect 
the FWHM(Ha) measurement, and, quite importantly, 
(3) the single-epoch mass estimates are based on H/3 line 
widths only. Therefore, the reverberation masses were re- 
computed using equation (1) based on FWHM(H 0), mea- 
sured in the mean and rms spectra, respectively, and the 
directly measured BLR sizes, Rblr, as listed in Tables 6 
and 7 by Kaspi et al. (2000). These H0 based reverber- 
ation masses are listed in Table 1 along with the single- 
epoch mass estimates. 

MBH(H/3,rms) is plotted in Figure 3a with the single- 
epoch mass estimates for the 18 objects in sample B (see 
Table 3; PG 1351+064 has no H(5 reverberation data). 
Similarly, A#bh(H/3, mean) is plotted in Figure 3b. The 
dashed line indicates a one to one relationship with no con- 
stant offset between the two mass measures. Some scat- 
ter is expected and observed and apart from the usual 
objects with deviating FWHM(H/3) measurements (the 
outlying objects are labeled), most of the scatter is likely 
caused by the intrinsic scatter in the BLR size — lumi- 
nosity relationship (see Figure 6 by Kaspi et al. 2000). 
This scatter is the cause of PG2130+099 being an out- 
lier in these diagrams. Disregarding PG1613+658, which 
all along has not behaved well, and for a moment the 
two objects with rather larger uncertainties in their mass 
estimates, PG1700+518 and PG1704+608, the scatter in 
the Mbh(H/3, rms)— Mbh(H/3, single-epoch) relationship 
is afterall comparable to the uncertainties in the mass es- 
timates. The best fit regression slopes are also plotted in 
Figure 3. The solid lines are the BCES bisectors when 
omitting the encircled object (PG1613+658). The dotted 
line is the BCES bisector to all the data points in each 
diagram. 

The results of the linear regressions and the detailed fit- 
ting parameters are listed in Table 4. The issue of estab- 
lishing the intrinsic relationship between these two mass 
estimates (single-epoch versus multi-epoch) is very similar 
to that of the most accurate representation of the intrinsic 
relationship between, say, the black hole mass, Mbh, and 
the stellar velocity dispersion in the host galaxy (e.g., Geb- 
hardt et al. 2000a; Ferrarese & Merritt 2000; Merritt & 
Fcrrarcsc 2001a). In both cases the relationship is sought 
between two variables with both measurement uncertain- 
ties and intrinsic scatter. Therefore, for the calibration 
of the single-epoch mass estimates, the use of the BCES 
linear regression algorithm is very important. 

The mass comparison of most importance for the cali- 
bration is that with Mbh(H/3, rms) (see § 2). Since both 
the single-epoch mass estimate and the reverberation mass 
are (different) measures of the same property, the BCES 
bisector is the most appropriate regression to use. It is 
apparent from Table 4, however, that using either of the 
BCES (Y|X) and bisector regressions yields the same basic 
results within the uncertainties. Namely, there is a pure 
unity relationship between Mbh(H/3, single-epoch), here- 
after M B h(H/3,S-E), and M BH (H/3,rms) within <0.25cr, 
for BCES(Y|X): 

log -Mbh (H/3 , rms) = (0.97 ± 0.13) x log M BB (R P , S - E) 
+ (0.05 ±1.08) (3) 
and for the BCES bisector: 
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log M B H(H/3,rms) 



(1.03 ± 0.14) x log M BH (H^, S - E) 
-(0.49 ±1.17) (4) 



The uncertainties used in the regression are, again, the 
symmetric errors on the logarithms based on the positive 
linear errors [i.e., = (loge) (<rx/X); § 4]. 

For completeness the single-epoch mass estimates are 
also compared with the reverberation masses based on the 
FWHM(H (3, mean) measured in the multi-epoch spectra 
(Fig. 3b). The outliers (labeled) are again the 'usual sus- 
pects' (§ 4.1.2 and above). Again, there is a unity relation- 
ship, but with less scatter (Table 4). This is expected as 
the single-epoch FWHM(H /?) also show less scatter with 
FWHM(H/3, mean) than with FWHM(H/3, rms) [§ 4.1]. 
When comparing the single-epoch masses to the masses 
computed by Kaspi ct al. (i.e., based on both Ha and 
H/3) larger scatter and poorer fits are obtained. This is 
expected, as explained, and stresses the need to use the 
reverberation masses based on H/3 measurements only. 

In conclusion, the calibration of single-epoch "optical" 
mass estimates to the multi-epoch reverberation masses is, 
with the current uncertainties, a one-to-one relationship 
with no significant zero-point offset and no apparent rea- 
sons to further correct Mbh(H/3, S-E). Note, however, that 
this assumes the line widths are measured as described in 
§ 4.1. The performance of the optical single-epoch mass 
calibration and how it compares to the UV single-epoch 
mass calibration, derived below, are discussed in § 7. 

6. CALIBRATION OF SINGLE-EPOCH MASS ESTIMATES: 
UV MEASUREMENTS 

In this section a calibration is determined of UV spec- 
tral measurements to reflect reasonable estimates of the 
central masses. Two subsets of data are used here. The 
primary sample (UVrev) is the collection of 26 AGNs and 
quasars with central mass determinations from reverber- 
ation mapping (Wandel et al. 1999; Kaspi et al. 2000) for 
which CiV line widths, FWHM(Civ), and UV continuum 
luminosities, ALa(1350A), are available in the literature. 
For this sample the UV measurements will be directly com- 
pared with the available reverberation mass determina- 
tions, Mbh(H/3, rms). By using this sample the UV cali- 
bration is based on fewer assumptions than if additional 
objects were included. This is because the single-epoch 
UV mass estimates are compared to independently mea- 
sured central masses and not to masses estimated based 
on, for example, optical measurements. 

Nevertheless, it is instructive to compare the results of 
the UVrev sample with those of a larger sample (named 
"sample UV" in the following). The latter sample con- 
sists of sample UVrev and 30 other PG quasars for which 
FWHM(C iv) and UV continuum luminosities are like- 
wise readily available (Table 2). For these additional 
quasars, no reverberation masses are available but the cen- 
tral masses are here determined from the calibration of 
single-epoch optical measurements, derived in § 5. It is 

5 Namely, only the measurements which deviated in excess of the assumed 10% measurement error were corrected. For quasars with strong 
and/or broad Fell emission but with no significantly strong, spiky narrow component, the BG92 measurement was considered the most accurate 
and fair representation of the FWHM(H/3) in the rms spectrum. Quasars which have 'uncorrected' FWHM(H/3) smaller than BG92 also have 
relatively strong narrow, sometimes spiky, components. In those cases the 'uncorrected' FWHM(H/3) was considered the most reasonable 
representation. 



clear that since these masses are estimates they are ex- 
pected to introduce some uncertainty, which is why they 
only serve to provide a comparison and an additional check 
on the performance of the calibrations. However, as dis- 
cussed later, this uncertainty appears to be within the 
scatter intrinsic to the reverberation masses, which indi- 
cates that the calibrations are relatively reliable. 

In § 4.1 it was necessary to remeasure the single-epoch 
FWHM(H j3) on the original, uncorrected data (from 
BG92) to ensure that the subtraction of the narrow emis- 
sion component by BG92 would not affect the line width 
measurement adversely. This was important for the use 
of FWHM(H/3) to estimate the central mass, as discussed 
there. Similarly, the single-epoch FWHM(H /3)'s were re- 
measured for the subset of 30 PG quasars, studied by 
BG92, which are used here to extend the primary sam- 
ple of reverberation mapped AGNs (UVrev). And sim- 
ilar to the approach adopted in § 4.1 each remeasurcd 
single-epoch FWHM(H/3) was compared to the BG92 
measurement and corrected based on the same consider- 
ations 5 . Similar to the UVrev objects the optical single- 
epoch continuum luminosities, ALa(5100A), were deter- 
mined from the tabulated spectral energy distributions of 
Neugebauer et al. (1987) or the photometry of Schmidt 
& Green (1983). The 'optical mass estimates' were then 
determined from equation (1) and (2), as justified in 
§ 5. The adopted optical measurements, FWHM(H/3) 
and ALa(5100A), and the single-epoch mass estimates, 
Mbh(H^), are listed in Table 1. 

The calibration of the UV measurements to reflect the 
central mass is described in the following. First, the mass 
is expected to depend on FWHM 2 (C iv). Anything else 
would not be physical (according to the virial theorem). 
Second, one can assume that the continuum luminosity 
enters to the power 0.7 as was established for the optical 
measurements (Kaspi et al. 2000; see also Appendix A for 
a discussion thereof), since the continuum luminosities at 
5100 A and 1350 A are generally expected to be related 
(via the assumed power- law relationship, L v oc u~ a ). The 
validity of this latter assumption can be tested by com- 
puting the BCES regressions between the already estab- 
lished masses and what can be called the normalized UV 
mass, nM BH ,uv= FWHM 2 (Civ) [A J L A (1350A)] 07 . The 
question is thus whether jiMbh.uv is related to Mbh to 
the first power and what is the global scaling factor? 
The regression results are presented in Table 5 for the 
two basic samples, namely sample UVrev and sample UV, 
and are discussed below. In this case where the rela- 
tionship between uJWbh.uv and the reverberation mass, 
Mbh(H/3, rms), will be used to calibrate other data it 
is very important to choose the most appropriate sta- 
tistical method to characterize the intrinsic relationship 
between these two variables. Different regressions [e.g., 
(Y|X), bisector, and orthogonal regressions] yields sta- 
tistically different results as they measure different mo- 
ments of the data (e.g., Isobe et al. 1990; Feigelson & 
Babu 1992). Which is then the most appropriate BCES 
regression (moment) to use for the calibration? This is 
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not clear even from a statistical point of view. There- 
fore, the following approach is adopted. First, BCES re- 
gression analyses are performed to test whether the rela- 
tionship between lognMeH,uv and log Mbh (H/3, rms) are 
indeed consistent with a slope of 1.0 to within the un- 
certainties. This justifies the assumption that Mbh.uv = 
constant x FWHM 2 (Civ) [AL A (1350A)]°- 7 as discussed 
above and also justifies the next step. Once this unity 
slope is established the problem reduces to the simple re- 
lationship with only one degree of freedom: y = x + a 
or a = y — x, where both y (= log Mbh(H/3, rms)) 
and x (= lognMeH.uv) have uncertainties. The constant 
zero point offset, a, can therefore be determined as the 
weighted mean of ' y — x ' determined for the indi- 
vidual objects. This approach limits the introduction of 
unnecessary systematic errors that would occur if both the 
slope and intercept are free parameters. Once a is deter- 
mined, it can be tested that the calibrated UV estimates, 
Mbh,uv= nMBH,uv x 10°, are indeed related to the estab- 
lished central masses, Mbh(H/3, rms), in a pure one-to-one 
relationship (i.e., a BCES bisector slope of unity and no 
zero-point offset). 

Figure 4 displays the distribution of ?tMbh,uv and 
Mbh (H/3, rms) for the UVrev sample (panel a) and for the 
full UV sample (panel b). The dotted line is the BCES 
bisector to all the displayed data points, while the dashed 
line represents the bisector when NGC4151 is excluded 
(see below). The solid lines are the BCES (Y|X) and 
(X|Y) regressions to the latter data. For the UVrev sample 
alone there is a larger difference between the bisector and 
BCES(Y|X) slopes (Fig. 4a; Table 5) than for the UV sam- 
ple (Fig. 4b), indicating the relatively larger scatter in the 
UVrev sample. [This is expected, however, as the 30 ad- 
ditional PC quasars have smaller intrinsic scatter in their 
optical mass estimates, Mbh(H/3), owing to their origin 
in the mass relationship in eqns. (1) and (2)]. As is also 
clear from these figures, NGC4151 has large uncertainties 
and is very much an outlier. Excluding this data point 
clearly has a significant effect on the slope in both cases 
(samples UVrev,b and UVb in Table 5). The BCES bisec- 
tor for sample UVrev, b then shows a slope of 1.0 to within 
0.3cr. The BCES (Y|X) regression is consistent with unity 
to within < 2.5cr; the uncertainty is somewhat large due to 
the intrinsic scatter in this sample. Note that, excluding 
also PG1704+608 and/or PG1613+658 does not signifi- 
cantly change the regressions. When the 30 additional PC 
quasars are included the statistical significance increases 
and also provides a larger mass range. In effect this allows 
the relationship to be better constrained (e.g., Fig. 4b). 
As seen in Table 5 for sample UV the (Y|X) slope is now 
increased to 0.7 and is improving further when NGC4151 
is excluded (sample UVb). In this case both the BCES 
(Y|X) and bisector slopes are consistent with a unity re- 
lationship to within at most ~ 2cr, which is acceptable. 
That is, both the UVrev, b and UVb samples show simi- 
lar relationships in the mass comparisons, although with 
different scatter and uncertainties. 

The fact that NGC4151 significantly changes the re- 
gression results argues that it should not be included in 
the calibration for the same reason that PG1613+658 was 
excluded from the optical calibration: it likely behaves dif- 
ferently than the bulk of the AGNs and quasars, and it is 



thus inappropriate to apply a calibration to large samples 
of AGNs that are based on and significantly affected by 
(possibly) strangely behaving objects. 

Since indeed log hMbh,uv and log Mbh (H/3, rms) are re- 
lated with slope = 1.0, the intercept can be determined as 
the weighted mean of the individual mass differences; their 
uncertainties are propagated from the individual errors on 
?t-Mbh,uv and Mbh (H/3, rms) . For sample UVrev, b the 
weighted mean is 

< log M BH (H/3, rms) -log nM BH ,uv > = 6.2±0.03 (5) 

where the uncertainty is the precision based on the propa- 
gated errors. The sample standard deviation is 0.45 show- 
ing the presence of real intrinsic scatter, as expected (e.g., 
Fig. 4). Excluding the other outliers (PG1613+658 and 
PG1704+608) does not significantly change the results. 
The weighted mean computed for the full UVb sample is 
the same, showing a difference only on the second signifi- 
cant digit (likewise for the precision and the sample stan- 
dard deviation). The fact that UVrev, b and UVb show 
very similar results confirms again that the optical single- 
epoch mass estimate is a fair representation of the rever- 
beration masses (based on the current data). 

How do the calibrated UV estimates, Mbh,uv= 
™Mbh,uv x 10 6 ' 2 , relate to the established central masses, 
Mbh(H/3, rms)? In Figure 5 Mbh,uv are plotted versus 
Mbh(H/3, rms) for the UVrev sample (panel a) and versus 
Mbh(H/3, rms) or M BH (H/3,S-E) for the full UV sample 
(panel b). The regression results (excluding NGC4151; 
Table 5) are also shown. The BCES (Y|X), (X|Y) (dot- 
ted lines) and bisector (solid line) regressions, and a pure 
1:1 relationship (dashed line) are also shown. The most 
important comparison here is that with the independently 
established reverberation masses in Figure 5a. The UV 
sample is shown in Figure 5b for completeness. The BCES 
bisector fit to the objects in the UVrev, b sample (excluding 
NGC4151) yields the following result: 

log M B H(H /3 ,rms) = 

(1.07 ± 0.14) x log M BH ,uv - 0.61 ± 1.11 (6) 

The bisector is consistent with a unity relationship to 
within 0.5(7, thereby indicating that the UV calibration 
is robust. 

When the additional 30 PG quasars are included 
(Fig. 5b) with single-epoch optical mass estimates a 
slightly steeper bisector is found. However, this sample is 
also consistent with a unity relationship in this calibration. 
In this case, it is so to within <~ 2a and is still acceptable. 
Although the two samples have different intrinsic scatter 
it is comforting to see that they yield consistent results, 
as pointed out earlier. 

In conclusion, the best calibration of the UV measure- 
ments appear to be 

log Mbh (Hp, rms or SE) = 

|"/ FWHM(CIV) \ 2 / A£ A (1350A) \° T 
° S 1000 km s- 1 J ^ 10 44 ergs s" 1 J 

+ 6.2 ±0.03 (±0.45) (7) 
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The last parenthesis contains the sample standard devia- 
tion of the weighted mean, which shows the intrinsic scat- 
ter in the sample. As opposed to the uncertainty or preci- 
sion of 0.03 on this weighted mean the standard deviation 
is probably more representative of the uncertainty in the 
offset for this (UVb) sample. It is also specific for the re- 
sult obtained in the logarithmic representation and makes 
little sense when linearized. Therefore, no uncertainties 
are listed for the linear scaling factor below. In linear rep- 
resentation, the calibration is: 



Mi 



BH,UV 



1.6 x 10' 



/ FWHM(CIV) \ 2 / A£ A (1350A) V 
^ 1000 km s" 1 J \10 i4 ergs s" 1 J 



(8) 



The mass calibrations rely strongly on the size — lu- 
minosity relationship and its intrinsic scatter dominate 
the uncertainties in the mass estimates. Such scatter is 
naturally expected given the nature of the objects and 
the nature of the obtained BLR sizes based on contin- 
uum variability (e.g., Netzer & Peterson 1997). Once a 
larger sample of AGNs with reverberation mapping masses 
is available, the -Rblr — L relationship should be updated. 
Because it matters how this relationship is established, 
a slightly modified approach to that adopted by Kaspi 
et al. (2000) is advocated in Appendix A. Note that, the 
modified approach does not significantly affect the current 
-Rblr — L relationship. 

7. HOW RELIABLE ARE THE CALIBRATIONS? 

It is of keen interest to ask how much of an error we typi- 
cally will make on the mass estimate using the virial mass- 
luminosity-FWHM relationship [eqn (1) and (2)] when in- 
stead of the FWHM(H /3,rms) and the mean monochro- 
matic luminosity, XL\, we use FWHM and \L\ measure- 
ments of a single-epoch spectrum, which is in fact a 'snap- 
shot' spectrum at any random given time? Clearly, this is 
only an approximation with the caveat that the estimated 
central mass may be off by a large factor. But what is the 
probability for that? In the following, the approximate un- 
certainties are briefly evaluated, assuming that the sample 
of 26 (18) nearby AGNs and quasars studied using rever- 
beration mapping techniques is representative for the UV 
(optical) calibration uncertainties. 

In Figure 6 the reverberation masses, Mbh (H/3, rms), 
are plotted against the deviation in the optical (panel a) 
and UV (panel b) mass estimates from this established 
mass. Offsets of 0.5 dex and 1.0 dcx from a perfect one- 
to-one relationship are indicated. The probabilities of es- 
timating the mass with a certain accuracy are summa- 
rized below and in Table 6. For the optical single-epoch 
mass estimates Figure 6a shows 17 objects out of the 18 
available to have single-epoch mass estimates deviating by 

6 The most well established reverberation mass is that of NGC 5548, which is based on many broad emission lines. Peterson & Wandel (2000) 
quote an uncertainty of ~0. 15 dex. In comparison, the propagated measurement errors on Mbh (H/3, rms) (Table 1) range from 0.04 dex to 1 dex 
(excluding those measurements with >100% errors) with an average propagated measurement error of 0.15 dex. The absolute uncertainty in 
the reverberation masses are difficult to access. If the measurement uncertainties of the stellar velocity dispersions arc assumed insignificant, 
the scatter in the M — <r relationship may indicate the typical uncertainty in the reverberation masses. There appears to be a general scatter of 
~0.15dex (Ferrarese & Merritt 2000), but it can be as much as 0.3 dex (Ferrarese et al. 2001). It is clear from Figure 2 by Ferrarese et al. that 
reverberation masses arc no more uncertain than the best masses derived from stellar kinematics. 



£1.0 dex. That is, there is a probability of ~95% of get- 
ting the mass accurate to within an order of magnitude 
using single-epoch optical spectrophotometry. Similarly, 
there is a ^90% probability of getting the mass accurate 
to within a factor of 6, and as much as a ~80% proba- 
bility of Mbh (H/3, S-E) being "correct" to within a factor 
of 3 (Table 6). In other words, the la uncertainty in 
Mbh (H/3, S-E) is a factor of <2.5. 

Three of the 26 objects in Figure 6b deviate in 
Mbh,uv(C IV, S-E) by more than 1 dex from the reverber- 
ation mass, Mbh(H/3, rms). There is, thus, a ^90% prob- 
ability of obtaining a central mass, accurate to within an 
order of magnitude, using single-epoch UV spectral mea- 
surements. As listed in Table 6 and illustrated in Fig- 
ure 6b, the probability of obtaining a mass to within a 
factor of 3 is as high as 70% (la error), and there is an 
even higher chance (85%) of the mass being within a factor 
of 6 based on the UV calibration presented here. 

The estimated uncertainty intrinsic to the reverberation 
technique is less than a factor 6 of 2. Therefore, it is ex- 
pected that the mass estimates based on the current opti- 
cal and UV calibrations are not all reliable within as little 
as a factor of 3. Yet it is reassuring that this relatively 
high accuracy can be achieved in a high 70% — 80% of the 
objects (UV and optical calibrations, respectively). 

Figure 6 shows a tendency of the lower-mass objects 
to be overestimated in the single-epoch masses for both 
the optical and UV calibrations. This may be an "arti- 
fact" of the objects with very narrow and spiky line pro- 
files. This is demonstrated in Figure 7 where the same 
mass deviations are plotted against FWHM(H f3, rms) - 
the width of the variable H /3 profile (left panel) and 
against FWHM(Civ, S-E) - the single-epoch Civ widths 
(right panel) — since a C iv rms width is only known for 
very few of the objects. Clearly, the most narrow-lined 
objects or those varying strongly in their narrow line core 
tend to be overestimated in their central masses. This is 
most pronounced for the optical mass estimates; note, it is 
inherit in the calibration process that the UV masses will 
scatter almost symmetrically about Mbh (H/3, rms). 

The good fortune is that the UV calibration appears to 
perform better for the higher mass and higher luminosity 
quasars, i.e., for FWHM(Civ) > 4000 kms" 1 and/or M B h 
> 5 x 10 7 M Q . This is very comforting as its primary appli- 
cation will be to samples of high-redshift quasars, which 
arc known to be of higher luminosity and are expected to 
occupy the higher end of the black-hole mass range ob- 
served (~ 4 x 10 7 - 1 x 10 9 M Q ; see also e.g. Laor 2000; 
Ho 2001). If the sample of PG quasars are representative 
most if not all quasars can thus have their central masses 
estimated to within a factor of ~5 or less (Figures 6 and 7). 

Since the mass estimates are very sensitive to the line 
width it is important that the signal-to-noise in the mea- 
sured spectra is sufficiently high that the FWHM is mea- 
sured relatively accurately (e.g., with a reliable uncertainty 
of 10% or less). This is worth keeping in mind, especially 
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when studying distant, faint quasars. It still applies, how- 
ever, that there is a 85% (90%) probability that the central 
masses can be estimated to within a factor of 6 based on 
UV (optical) data. This is reliable enough for many sta- 
tistical applications. 

In conclusion, both the optical and UV single-epoch 
mass estimates indicate approximately similar probabili- 
ties of obtaining reliable mass estimates with the optical 
calibration performing slightly better in general. The 1 a 
uncertainty is a factor of 3 or better. The UV calibration 
performs very well for the higher mass, broader lined ob- 
jects where practically all the mass estimates are within a 
factor of ^5 or less of the reverberation mapping masses. 

8. SUMMARY AND CONCLUSIONS 

Virial estimates of central black hole masses based on ei- 
ther optical or UV single-epoch spectrophotometry are cal- 
ibrated to recently measured masses of nearby AGNs and 
quasars using reverberation mapping techniques. Single- 
epoch spectral measurements allow the central masses of 
distant AGNs and quasars to be easily estimated. This 
is important when more direct mass measurement tech- 
niques cannot easily be applied; moreover, this method 
allows estimates of masses of large samples of AGNs in a 
short time span. The following conclusions are reached: 

• The current data show that the above-mentioned 
mass estimates are best made as follows: 

— The signal-to-noise in the single-epoch spec- 
tra is high enough that the FWHM can be 
measured reliably to an accuracy of 10% or 
better. 

— The FWHM is measured on a H/3 or Civ 
emission line profile corrected for strong, con- 
taminating Fen emission. But, for higher lu- 
minosity sources (i.e., quasars), narrow- line 
subtraction should not be attempted. 

— The continuum luminosity measured at 
5100A or 1350A is used to estimate the rele- 
vant distance of the H (3 or C iv emission line 
gas, respectively, using the empirical -Rblr — 
L relationship (Kaspi et al. 2000). 



>4000 kins" 1 essentially all masses here are ac- 
curate to within a factor of 3 to ^5. Thus these 
calibrations seem to perform very well where they 
are most needed, that is at high redshift for higher 
luminosity, more massive, quasars. 

• The currently obtainable accuracy and associated 
probabilities are quite fair given the intrinsic un- 
certainty (factor less than 2) in the reverberation 
masses. An increase in the accuracy of the cali- 
brations require a smaller scatter in the BLR size 
— luminosity relationship. A decrease in the mea- 
surement uncertainties of the BLR size from rever- 
beration studies will help. However, it is equally 
important that the BLR size be determined for as 
large a sample of AGNs and quasars as possible 
(and larger than the current sample of 34 objects) 
spanning a range of AGN properties. 
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• The 1 a uncertainty in the optical and UV calibra- 
tions are a factor of 2.5 and 3, respectively. 

— Mass estimates based on single-epoch optical 
spectra are statistically consistent with the 
reverberation masses to within the uncertain- 
ties for the current database. The 'optical 
single-epoch mass estimate' measure the re- 
verberation masses to within factors of 3, 6, 
and 10 with probabilities of 80%, 90%, and 
95%, respectively. 

— The 'UV single-epoch mass estimate' measure 
the reverberation masses to within factors of 
3, 6, and 10 with probabilities of 70%, 85%, 
and 90%, respectively. 

• The most deviating mass estimates are found for 
the lower mass, more narrow- lined and/or lower lu- 
minosity AGNs. For quasars with FWHM(Civ) 
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Fig. 1. — The distribution of FWHM(H/3) in the mean 
(panel a) and rms (panel b) multi-epoch spectra and single-epoch 
FWHM(H P) of the PG quasars presented by Kaspi et al. (2000) 
and MrkllO and Mrk335 (Wandel et al. 1999). The dotted line 
indicates a pure one-to-one relationship. The objects discussed in 
the text are labeled. The open squares denote Seyfert Is while tri- 
angles show measurements for the quasars. The three solid circles 
show the measurements of Boroson & Green (1992) which arc based 
on spectra with both Fe II emission and the narrow core component 
subtracted. These FWHM measurements arc not always represen- 
tative of the BLR velocity dispersion needed for this study (see 
text). The solid line is the best fit BCES bisector regression line 
based on all the objects in the diagram except PG1613+658 (sam- 
ple C; Tabic 3). The BCES (Y|X) and (X|Y) regressions are not 
plotted as they crowd the bisector. The single-epoch FWHM(H/3) 
scatter around a one-to-one relationship with FWHM(H /3, mean) 
to within 15% — 20% variation and around a similar relationship 
with FWHM(H/3, rms) to within 20% - 25%. Note, the ordinate 
range is different in the two diagrams. 

Fig. 2. — (a) The distribution of the mean AL A (5100A) multi- 
epoch measurements (Wandel et al. 1999; Kaspi et al. 2000) based 
on monitoring data with respect to the single-epoch AL.\(5100A) 
measurements of Neugebauer et al. (1987) and Schmidt & Green 
(1983). The errors in the mean AL A (5100A) arc the rms around 
this mean. The errors in the single-epoch AL A (5100A) are propa- 
gated errors (see text). The short-dashed line (centrally positioned) 
denotes a one-to-one relationship. The dotted and long-dashed lines 
represent ±30% and ±60% luminosity variations, respectively (see 
text), (b) The best fit regression lines are shown for the BCES 
bisector (solid line). The BCES (Y|X) and (X|Y) regression lines 
would crowd the bisector, if plotted. All these BCES regression fits 
are consistent with a slope of 1.0. The single-epoch luminosities are 
offset by +0.138 dex at AL A (5100A) «44.8 ergs s _1 , the mid-range 
luminosity, based on the BCES bisector. 

Fig. 3. — The reverberation masses derived from the rms (panel 
a) and mean (panel b) spectra plotted versus the single-epoch mass 
estimates based on optical spectral measurements. Triangles denote 
quasars, while squares denote Seyfert Is. The dashed line indicates 
a unity relationship. The dotted line is a BCES bisector regression 
line to all the PG quasars with H (i measurements (shown), while the 
solid line is the BCES bisector when PG1613 is excluded (see text). 
As expected, Mbh(H/3, S-E) show less scatter with Mbh(H0, mean). 
Both relationships are consistent with a one-to-one relationship to 
within the errors. 



Fig. 4. — Established and estimated central masses based on 
optical data plotted versus the UV measurements. This distribu- 
tion is the basis of the calibration of the UV measurements. The 
"optical masses" (ordinate) in panel (a) consist of the reverberation 
masses, Mbh(H/3, rms), (filled, circled squares) derived from the 
rms spectrum and based on H (i only. In panel (b) these masses are 
supplemented with the single-epoch mass estimates, Mbh(H/3, S-E), 
(filled squares) for the 30 PG quasars with no reverberation map- 
ping masses. Dotted lines: BCES bisector regression lines to all the 
objects in each diagram. Solid lines: BCES(Y|X) and (X|Y) regres- 
sions to all the objects except NGC4151 (see text). Dashed lines: 
the BCES bisector for all objects excluding NGC4151. 

Fig. 5. — The central mass estimates based on the calibra- 
tion of the UV spectral measurements are compared to central 
masses measured (and/or estimated; panel b) based on optical data. 
Dashed line: pure unity relationship. Dotted lines: BCES(Y|X) 
and BCES(X|Y) regression lines. Solid line: the BCES bisector. 
NGC4151 was excluded from the regression analysis, (a) Mbh,uv 
estimates versus the established central masses, M^n(H(3, vms) 
[UVrcv sample only], (b) Mbh.uv estimates arc plotted for the 
full UV sample versus the "optical masses" described in Figure 4b. 
For both diagrams the mass relationships arc consistent with a unity 
relationship within the uncertainties. 

Fig. 6. — The established central masses, Mbh(H/3, rms), based 
on optical multi-epoch spectral measurements plotted versus the 
deviations of mass estimates based on calibrated single-epoch spec- 
tra, (a) Deviations of masses estimated from optical spectra, i.e., 
M B h(H/3, S-E) divided by M B h(H/3, rms). (b) Deviations in the 
masses based on UV spectral measurements, Mbh,uv(C iv). The 
uncertainties in the abscissa are the (propagated) uncertainties in 
the single-epoch masses (i.e., not the mass deviation error). A 
strictly unity relationship is indicated by the solid line. Offsets of 
±0.5 dex (±ldex) are indicated by the dotted (dashed) lines. 

Fig. 7. — The mass deviations from Figure 6 plotted here versus 
the H (i and C IV line widths. The optical mass deviations are plot- 
ted against FWHM(H (3, rms), the width of the rms profile (i.e., the 
variable part), in (a) and versus the single-epoch FWHM (Civ) in 
(b). Lines and symbols are as in Figure 6. Notice that the larger 
mass discrepancies tend to occur for the most narrow lined objects 
in both cases (or those with the strongest variability occurring in 
the narrow line core: PG1613±658 and PG1704±608). 



APPENDIX 

A. THE BLR SIZE — LUMINOSITY RELATIONSHIP REVISITED 

The aim of this section is twofold, namely (1) to demonstrate the importance of using regression analysis that is 
appropriate for the nature of the data set, and (2) to advocate a careful assessment of the object sample on which 
important calibrations, such as the -Rblr — L relationship, are based. While this approach yields an Rblr — L relationship 
consistent with equation (2) for the current object sample and data base, it is expected to become important once more 
data are available, allowing better constraints to be placed on this relationship. 

Kaspi et al. (2000) established the empirical relationship [equation (2)] between the BLR size as measured by rever- 
beration mapping and the intrinsic continuum luminosity of the objects using linear regression techniques that take 
errors in both variables into account. However, as the relationship is sought between two variables that have both mea- 
surement uncertainties and intrinsic scatter, it is important to use the bivariate correlated errors and intrinsic scatter 
(BCES) algorithm (Akritas & Bcrshady 1996) to establish the best size — luminosity relationship. The results of apply- 
ing this more appropriate linear regression analysis is shown in Table Al. The regression was done between log -Rblr. 
and log[ALA(5100A)/10 44 ergs s _1 ] to allow a linear regression. Although i?BLR is physically speaking a function of the 
luminosity and thus the BCES(Y|X) is the appropriate regression to use, the BCES(X|Y) and bisector regression results 
are also listed for completeness. The more accurate BLR size of 5.9^2'q f° r NGC4051 (Peterson et al. 2000) is used here 
instead of that listed by Kaspi et al. (2000), but is not the cause of the differences discussed below. 

When intrinsic scatter is accounted for the luminosity clearly enters with a power different from 0.7 (Table Al), the 
luminosity power established by Kaspi et al. The linear equivalents of the log — log BCES regressions are for the sample 
including NGC 4051: 



i?BLR = (20.9 ±4.7) 



AL A (5100A) 



-, (0.38±0.14) 



10 44 ergs s- 1 



(Al) 
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for the Seyfert Is alone, and for the PG quasars alone: 



Rblr = (41.8 ±8.3) 
The linear relationship for the combined sample is: 

i? BL R = (32.4 ± 4.7) 



AL A (5100A) 



10 44 ergs s" 



AL A (5100A) 



10 44 ergs s- 1 



(0.52±0.09) 



(0.58±0.09) 



(A2) 



(A3) 



The uncertainties used in the regression are the symmetric errors on the logarithms based on the positive linear errors. 
The errors on the proportionality factor are the linearized errors of the errors on log -Rblr ; listed for comparison with 
eqn. 2. 

When the intrinsic scatter in the data is accounted for the luminosity dependence of the BLR size weakens and the 
uncertainties are increased [eqn. (A3)]. However, these larger uncertainties are perhaps more representative than those 
returned by a weighted regression alone (eqn. 2). Note that the slopes of the individual samples [eqn. (Al) and (A2)] are 
consistent within their errors. 

Peterson ct al. (2000) find that NGC 4051 is peculiar as it has a larger BLR size for its luminosity than the other objects, 
deviating by almost 3a. Excluding this object they find a luminosity dependence to the power 0.6 using a variance weighted 
regression algorithm similar to Kaspi et al., who included all objects in their regression. When NGC 4051 is excluded the 
BCES linear regressions change as follows: 



for the sample of Seyferts 



Rblr = (23.4 ± 4.5) 



-Rblr = (30.2 ± 5.6) 



AL A (5100A)l (0 - 56± °- 26) 



10 44 ergs s" 



AL A (5100A) 



10 44 ergs s" 



(0.66±0.09) 



(A4) 



(A5) 



and for the combined PG and Seyfert sample. 

The value for NGC 4051 is seen to have a significant effect on the slope. Excluding this object makes the Seyfert 
sample regression consistent with that of the combined sample to within the larger uncertainties. It is also consistent 
with the result of Peterson et al. (2000). Both the modified regression slopes, shown in their linear representations in 
equations (A4) and (A5), also agree with the result by Kaspi et al. (2000) within the uncertainties. Which slope should 
be used to estimate the central masses? Similar to the mass calibrations (§§ 5 and 6) one can argue that a calibration like 
the present should not include an extreme outlier that may not be representative of the common AGN, especially when it 
has a significant effect on the result. Therefore equation (A5) seems the most appropriate size — luminosity relationship 
to use. However, given the current scatter and uncertainties, using a slope of 0.7 and a scaling factor of 32.9 makes little 
difference to the mass estimates. However, the larger uncertainties in equation (A5) are more representative than those 
of eqn. (2). 

Note, the regression results do not change significantly from those listed when the linear negative errors on the BLR 
size is used instead of the positive ones, or if the largest of the two is used for each object; this in fact illustrates the 
robustness of the BCES regression. 

The precise value of the luminosity power in the size — luminosity relationship is important for the single-epoch mass 
calibrations. Mass estimates of distant AGNs are highly dependent thereon as the masses cannot at present be determined 
independently in a direct manner. This stresses the dire need to obtain BLR size measurements for a larger sample of 
AGNs and to a higher degree of accuracy such that this relationship can be further constrained and studied. Once 
such data are available the approach presented here (use of BCES regression analysis and exclusion of objects, such as 
NGC 4051, which are not good representations of the typical AGN population) should be adopted. 
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Table Al 

BCES Regression for BLR Size - Luminosity Relationship 



Independent 


Dependent 


Sam- 


N 


Slope 


Intercept 


Variable 


Variable 


ple 


(#) 


ierror 


ierror 


Bivariate Errors and Intrinsic Scatter Regression (BCES) 


logR a 


logAL A b 


All 


34 


0.58±0.09 


1.51±0.06 


log XL\ 


logR 


All 


34 


0.63±0.11 


1.51±0.06 


BCES bisector 




All 


34 


0.61±0.09 


1.51±0.06 


logR 


log AL A 


Sy Is 


17 


0.38±0.14 


1.32±0.10 


log AL A 


logR 


Sy Is 


17 


0.56±0.25 


1.41±0.15 


BCES bisector 




Sy Is 


17 


0.47±0.17 


1.36±0.11 


logR 


log AL A 


PCs 


17 


0.52±0.09 


1.62±0.09 


log AL A 


logR 


PGs 


17 


0.45±0.32 


1.67±0.21 


BCES bisector 




PGs 


17 


0.49±0.20 


1.64±0.18 


logR a 


logAL A b 


A' 


33 


0.66±0.09 


1.48±0.07 


log AL A 


logR 


A' 


33 


0.71±0.12 


1.47±0.06 


BCES bisector 




A' 


33 


0.68±0.10 


1.47±0.06 


logR 


log AL A 


S' 


16 


0.56±0.26 


1.37±0.10 


log AL A 


logR 


S' 


16 


0.77±0.26 


1.46±0.15 


BCES bisector 




S' 


16 


0.66±0.19 


1.41±0.10 



Sample All: All Seyfert Is and PG quasars (ref [1,2]) 
Sample Sy Is: The Seyfert Is only (ref [1,2]) 
Sample PGs: The PG quasars (only) (ref [2]) 

Sample A': All Seyfert Is and PG quasars except NGC4051. See text 
Sample S': The Seyfert Is except NGC4051. See text 



a log[(BLR Size, Rblr)/ light-days] 
b log[AL A (5100A)/ 10 44 ergs s" 1 ] 

References. — (1) Wandel ct al. 1999; (2) Kaspi et al. 2000 
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Table 1 

Optical Spectral Parameters and Masses 



Object 



Alt. 

Name 



Red- 
shift 



FWHM 

(H/3,S-E) 



log[AL A ] b 



log [M/Mq] c 
(H/3.S-E) 



log [M/M ] a 
(H/3, rms) 



log [M/Mq] c1 Ref. 
(Hj3, mean) 



3C 120 
3C 390.3 
Akn 120 
Fairall 9 
Mrk79 
Mrk509 
NGC 4151 
NGC 3783 
NGC 5548 
NGC 7469 



0.033 
0.057 
0.033 
0.046 
0.022 
0.035 
0.003 
0.010 
0.017 
0.017 



7 48 + ' 21 

-0.27 

8 57 + ' 12 
8 j7 -0.21 
„ +0.08 
s ' 27 -0.12 

7 92 + ' 11 
-0.32 

8.01 tolt 

7 96 + 05 
' aD -0.06 
- ns +0.23 
7 ' 08 -0.38 
+ 0.30 
-0.96 
+ 0.07 
-0.07 
+ 0.30 
6.88 



7.04 
7.97 
6.88 



+ 0.22 
-0.28 
+ 0.12 
-0.21 
+ 0.08 
-0.12 

7 90 + ' 11 
, au -0.31 
+ 0.14 
-0.34 
+ 0.05 
-0.05 
+ 0.23 
-0.38 
+ 0.30 
-0.97 
+ 0.07 
-0.07 
R „, +0.30 
6 ' 81 6.81 



7.36 
8.53 
8.26 



7.72 



6.97 
8.09 



PG 0003+199 
PG 0007+106 
PG 0026+129 
PG 0050+124 
PG 0052 + 251 
PG 0157+001 
PG 0804+761 
PG 0844+349 
PG 0921 + 525 
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PG 1259+593 
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PG 1351 + 640 
PG 1352 + 183 
PG 1411 + 442 
PG 1415 + 451 
PG 1416-129 
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PG 1444+406 
PG 1501 + 106 
PG 1512 + 370 
PG 1534+580 
PG 1545 + 210 
PG 1613 + 658 
PG 1704+608 
PG 2130 + 099 
PG 2209 + 184 
PG 2251 + 113 
PG 2308 + 098 
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III Zw 2 
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Mrk 1014 
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Mrk 110 
Mrk 705 
K348-7 
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PKS 
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PKS 
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Mrk 1383 
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4C 37.43 
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3C 323.1 
Mrk 876 
3C 351 
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PKS 
4C 09.72 
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0.089 
0.142 
0.061 
0.155 
0.164 
0.100 
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0.029 
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0.425 
0.144 
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0.114 
0.129 
0.086 
0.077 
0.267 
0.036 
0.371 
0.030 
0.266 
0.129 
0.371 
0.061 
0.070 
0.323 
0.432 



1640 


± 


164 




44 


00 


± 





13 


5 


7 


14 


± 





12 


6.58 


+ 0.14 

-0.14 


6.80 


+ (> 

-0 


12 

09 


1 


3979 


± 


398 




44 


59 


± 





07 


4 


8 


29 


± 





10 




+ 0.08 
-0.11 




+ 
-0 


08 

11 




1860 


± 


250 


2 


44 


90 


± 





07 


4 


7 


85 


± 





12 


7.48 


7.86 


6 


1240 


± 


124 


2 


44 


52 


± 





12 


5 


7 


26 


± 





11 




+ 0.10 
-0.10 




+0 
-0 


10 
10 




5200 


± 


520 


2 


44 


87 


± 





07 


4 


8 


72 


± 





10 


8.61 


8.49 


6 


2460 


± 


320 


2 


44 


81 


± 





07 


4 


8 


03 


± 





11 




+ 0.04 
-0.04 




+ 
-0 


04 
04 




3070 


± 


307 


2 


44 


93 


± 





07 


4 


8 


31 


± 





10 


8.13 


8.31 


6 


2420 


± 


242 


2 


44 


34 


± 





07 


4 


7 


69 


± 





10 


7.45 


+ 0.15 
-0.21 


7.42 


+ 
-0 


15 
21 


6 


1816 


± 


182 


3 


43 


55 


= 





12 


5 


6 


91 


± 





12 


6.89 


+ 0.14 
-0.21 


6.75 


+ 
-0 


13 
19 


1 


1742 


± 


248 


3 


43 


66 


= 





13 


5 


6 


95 


± 





14 




+ 0.06 
-0.09 




+ 
-0 


06 
09 




3130 


± 


470 


2 


45 


30 


± 





06 


4 


8 


58 


± 





13 


8.22 


8.27 


6 


2640 


± 


264 


2 


44 


83 


± 





07 


4 


8 


11 


± 





10 














5360 


± 


536 


2 


45 


43 


= 





06 


4 


9 


14 


± 





10 














4973 


± 


497 


3 


45 


41 


= 





06 


4 


9 


07 


± 





10 














6190 


± 


619 


2 


45 


47 


= 





06 


4 


9 


29 


± 





10 














4570 


± 


457 


2 


44 


59 


= 





07 


4 


8 


41 


± 





09 














2920 


± 


292 


2 


45 


27 


= 





07 


4 


8 


50 


± 





10 














1820 


± 


182 


2 


44 


00 


= 





08 


4 


7 


20 


± 





10 














3715 


± 


372 


3 


44 


42 


= 





07 


4 


8 


12 


± 





09 




+ 0.10 
-0.16 




+ 
-0 


09 
15 




1860 


± 


186 


2 


44 


94 


= 





07 


4 


7 


88 


± 





10 


7.51 


7.70 


6 


5190 


± 


1020 


2 


45 


50 


= 





06 


4 


9 


17 


± 





15 




+ 0.06 
-0.06 




+ 
-0 


06 
06 




3520 


± 


352 


2 


45 


89 


= 





07 


4 


9 


10 


± 





11 


8.63 


8.82 


6 


3360 


± 


336 


2 


44 


27 


= 





07 


4 


7 


93 


± 





10 


7.95 


+ 0.17 
-0.27 


7.94 


+ 
-0 


17 
27 


6 


830 ± 


83 


2 


43 


67 


= 





08 


4 


6 


29 


± 





10 














3390 


± 


1139 


2 


45 


77 


= 





06 


4 


8 


98 


± 





23 














3400 


± 


340 


2 


45 


70 


= 





06 


4 


8 


94 


± 





10 




+ 0.14 
-0.46 




+ 
-0 


14 
46 




5320 


± 


532 


2 


44 


86 


= 





07 


4 


8 


73 


± 





10 


8.70 


8.50 


6 


2940 


± 


476 


2 


44 


84 


= 





07 


4 


8 


21 


± 





13 














5660 


± 


566 


2 


44 


71 


± 





07 


4 


8 


69 


± 





10 














3600 


± 


360 


2 


44 


68 


= 





07 


4 


8 


27 


± 





09 




+ 0.14 
-0.20 




+ 
-0 


14 
20 




2670 


± 


267 


2 


44 


50 


= 





07 


4 


7 


89 


± 





10 


8.05 


7.96 


6 


2620 


± 


262 


2 


44 


42 


= 





07 


4 


7 


81 


± 





09 














3766 


± 


377 


3 


44 


95 


= 





07 


4 


8 


50 


± 





10 














5940 


± 


594 


3 


44 


74 


= 





07 


4 


8 


75 


± 





10 


8.63 


+ 0.13 
-0.24 


8-73 


+ 
-0 


13 
24 


6 


1450 


± 


145 


2 


44 


42 


= 





07 


4 


7 


30 


± 





10 














2480 


± 


248 


2 


45 


08 


= 





06 


4 


8 


23 


± 





10 














3465 


± 


347 


3 


44 


17 


= 





09 


4 


7 


88 


± 





10 














4352 


± 


435 


3 


45 


41 


= 





06 


4 


8 


95 


± 





10 














3060 


± 


306 


3 


43 


58 


= 





09 


4 


7 


36 


± 





10 














5961 


± 


596 


3 


45 


26 


= 





06 


4 


9 


12 


± 





10 




+ 0.18 
-0.20 




+ 
-0 


18 

20 




8450 


± 


845 


2 


44 


65 


= 





07 


4 


8 


99 


± 





10 


7.55 


8.45 


6 


1180 


± 


150 


3 


45 


50 


= 





06 


4 


7 


88 


± 





12 


6.87 


+ 0.26 
-1.12 


7.57 


+ 

- 1 


27 
03 


6 


2330 


± 


233 


2 


44 


37 


= 





07 


4 


7 


68 


± 





10 


8.42 


+ 0.13 
-0.07 


8-23 


+ 
-0 


13 
07 


6 


6500 


± 


650 


2 


44 


30 


= 





07 


4 


8 


52 


± 





10 














4160 


± 


710 


2 


45 


48 


= 





06 


4 


8 


96 


± 





14 














7920 


± 


792 


2 


45 


52 


= 





11 


5 


9 


57 


± 





12 















a FWHM(Hfj) measured in the single-epoch spect 
b log[AL A (5100A)/crg-s s" 1 ]. 
c The central mass (and uncertainties) est: 



"The cen 
nd R BLR : 



al : 



.ss (a 



s) dctc 
rcstfrai 



atcd based 
•mined fron 



m single-epoch optical spect: 



nulti-epoch spectrophotometry: FWHM(H/3) 
ure me lit;- arc determined u;diifr reverberation 



spectrum, respectively, 



References. — (1) Wandel, Peterson, & Malkan 1999; (2) Boroson &; Green 1992; (3) This author's own measurements of the original, uncorrected 
spectra presented by ref.[2]; (4) Neugebaucr ct al. 1987; (5) Schmidt & Green 1983; (6) Kaspi ct al. 2000. 
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Table 2 

Ultraviolet Spectral Parameters and Masses 



Object 


Redshift 


FWHM 1 


Ref. b 


log [AL; 




log [M/M Q ] d 






(CIV, S-E) 










(CIV, 


S-E) 


IP i on 
o 1ZU 


n 033 

U.UOO 


QQ07 -t- ^so 

OOU ( ZE OOU 


1 

1 


44.17 


± 0.09 


7.48 


± 0.12 


o^ oyu.o 


0^7 
U.UO 1 


7700 -t- 7^0 
/ / UU ZH 1 ou 


9 


43.74 


± 0.17 


7.79 


± 0.17 


Akn 120 


0.033 


4908 ± 802 


1 3 


44.70 


± 0.12 


8.07 


± 0.17 


lairall 9 


0.046 


3826 ± 805 


1,3 


44.88 


+ 


0.09 

1) 1 1 


7.98 


± 0.20 


Mrk 79 


0.022 


5552 ± 560 


1 


43.87 


± 0.11 


7.60 


± 0.13 


Mrk 509 


0.035 


5303 ± 658 


1,3 


44.46 


+ 


0.08 
0.09 


7.97 


±0.13 


NGC4151 


0.003 


2190 ± 562 


3 


41.36 


+ 


0.29 

n.fin 


5.03 


+ 


0.37 
n .32 


NGC3783 


0.010 


3476 ± 336 


4 


43.10 


± 0.15 


6.65 


± 0.15 


NGC 5548 (mean) 


0.017 


6422 ± 2529 


1,3,5,6 


4o.D4 


+ 0.11 
-0. 16 


7.56 


± 0.36 


— TUR ( 1 9891 


0.017 


5520 ± 380 


5 


43.51 


± (0.10) 


7.34 


± 0.10 


_ ffQ r T (1 QQQ^ 

— .no i i lyso j 


0.017 


8950 ± 570 


6 


43.43 


± (0.10) 


7.70 


± 0.10 


NGC 74fi9 


0.017 


4570 ± 380 


7 


43.49 


± 0.08 


7.16 


± 0.10 


pr; 000^4-1 QQ 


0.025 


3777 


± 


37^ 
O f 




44.23 


4_ 


10 


7.52 


_l_ 

□z 


19 
U. 1Z 


pp ono7-ui nfi 

I V_j UUU I TlUU 


08Q 

u.uoy 


4390 


± 


730 
( OU 


Fl 8 LI 


44.72 


4- 


11 


7.99 


-L- 
-E 


17 
U. 1 1 


PP 009R4-1 9Q 
i v_j uuzu^izy 


142 


3261 


± 


^97 

OZ / 


PI Q T 1 

r i ,o,-Lj1 


45.11 


-(- 


07 


8.00 




n it: 

U. lO 


pp oo^04-i 94 


0fi1 


1885 


± 


1 QO 

iyu 


Q 
y 


44.65 


4_ 


0^ 
U.UO 


7.21 


-L- 

m 


10 
U. 1U 


PP 00^94-9^1 


U. -LOO 


6653 


± 


uuo 




45.29 


4- 

4Z 


12 


8.75 


-E 


(1 I 1 ) 

U. lo 


PP 01 ^74-001 

I V. J UIO 1 tUWI 


1 fi4 

U. -LU4 


4184 


± 


420 




45.07 


4- 


OQ 
u.uy 


8.19 


4- 

-E 


12 


PP 0R044-7K1 


1 00 
U. -LUU 


4241 


± 


90*3 


ri a t i 

r i ,o,-Lj1 


45.07 


4_ 

ZL 


OS 
U.UO 


8.20 


4: 


OS 
U.UO 


PP 0R44-U34Q 


0fi4 

U.UU4 


4439 


± 


1 4Q 

i^y 


ri a t i 

r i ,o,-Lj1 


44.57 


4_ 
ztz 


01 
U.Ul 


7.89 


4- 
4: 


0^ 
U.UO 


PP 0Q91 4-^9^ 


03^ 
U.UOO 


3782 


± 


OOU 




43.71 


4_ 

ZL 


1 fi 
U. ID 


7.15 


4- 
4: 


17 
U. 1 / 


PP 009^4-1 9Q 


09Q 


4103 


± 


410 




43.74 


4^ 


0.15 


7.24 


j- 


0.15 


PP 00^4-41 4 
i v_j uyooT'ii'i 


93Q 

u.zoy 


3067 


± 


^i n 


1 

1U 


45.51 


4_ 
ztz 


(U.10J 


8.23 


4- 
4: 


19 
U. 1Z 


PP 1 01 94-008 


0.185 


5816 


± 


120 


Fl 8 LI 


44.82 


4- 

ZIZ 


OQ 
u.uy 


8.30 


4- 

EE 


OS 
U.UO 


PP 1 04Q 00^ 
x vj lufiy — uuo 


3^7 
U.OO I 


3675 


± 


4^0 

4:OU 


Fll L14 


45.75 


4- 

ZIZ 


1 1 

U. J. 1 


8.56 


4- 

EE 


14 


PP 1 1 004-779 
r v_j iivjut^ i f z 


313 


8664 


± 


QOO 

yuu 


PI Q T 1 

r i ,o,-Lj1 


45.60 


4_ 

ZL 


OQ 

u.uy 


9.20 


4_ 

EE 


19 


PP 1 1 1 } OOfi 
1 V_j _L_LUO — UUU 


49^ 

U.4ZO 


3500 


± 


4^0 


Fll L14 


45.87 


4_ 

ZL 


Oil 
U. 11 


8.60 


4_ 

EE 


n it; 
U. 10 


PP 1 1 144-44^ 


0.144 


2246 


± 


225 


1 


44.70 


4- 
ztz 





7.39 


4- 

EE 


^7 
U.O 1 


PP 111 fi-l-91 ^ 
1 V? _L 1 IUtZ _LO 


1 77 
U. 1 it 


5329 


± 


ooo 


PS 1 9 T 1 9 

-T 0,1Z,Jj1Z 


45.61 


4- 

ZL 


f 1 0~1 
(U.1UJ 


8.78 


4- 

EE 


U. 1U 


PP 111 Q_l_ 1 90 
i v_j iiiy^izu 


04Q 

u.u^y 


5425 


± 


040 


I 


44.09 


4_ 

ZL 


10 
u. 1U 


7.73 


4- 
4: 


19 


PP 1 9094-9R1 


U. -LUO 


2941 


± 


9Q0 
zyu 


PS 1 9 T 1 9 


44.27 


4- 

ZIZ 


^0 1 0~1 

(^u.iu; 


7.33 


4- 

EE 


12 


PP 1 91 1 4-1 43 


08^ 
U.UOO 


2981 


± 


247 


PI O T 1 

r i ,o,-Lj1 


45.01 


± 


0.08 


7.86 


± 


0.10 


PP 1 91 R4-0RQ 


334 

U.004 


3039 


± 


300 


12 


45.62 


± 


(0.10) 


8.30 


± 


0.12 


PP 1 99fi4-093 


U. -LOO 


4141 


± 


381 


PR 1 T 1 

r o, ±u,±j-Lu 


46.34 


± 


(0.10) 


9.07 


± 


0.11 


PP 1 99Q4-904 


0fi4 

U.UU4 


4904 


± 


490 




44.59 


± 


0.08 


7.99 


± 


0.11 


PP 1 944-U09R 


048 
U.U40 


2531 


± 


255 




43.69 


± 


0.13 


6.79 


± 


0.14 


PP 1 9^Q4_^Q'? 


479 


5925 


± 


800 


PI 1 T 1 4 


46.05 


± 


0.11 


9.18 


± 


0.15 


PP 1 309 — 1 09 


0.286 


3634 


± 


365 


12 


45.90 


± 


(0.10) 


8.65 


± 


0.12 


PP 1 3074-08^ 
1 V i _LOU I tUOO 


U. -LOO 


6399 ± 


1854 


PI Q T 1 

r i ,o,-Lj1 


45.12 


± 


0.10 


8.60 


± 


0.26 


PP 1 ^OQ-I-^^^ 


0.184 


4962 


± 


500 




44.87 


± 


0.34 


8.20 


± 


0.37 


PP 1 3^1 4-fi40 


087 
U.UO I 


2864 


± 


751 


PI Q T 1 

r i ,o,-Lj1 


44.65 


± 


0.06 


7.57 


± 


0.23 


PP 1 3^94-1 83 


U. -LOO 


3927 


± 


395 




44.91 


± 


0.09 


8.03 


± 


0.12 


PP 1 41 1 4-449 
I V_j L^tL J.-J-44Z 


08Q 

u.uoy 


2617 


± 


260 




44.46 


± 


0.08 


7.36 


± 


0.11 


PP 1 41 ^4-4^1 


114 


6768 


± 


680 




44.36 


± 


0.13 


8.11 


± 


0.14 


PP 1 41 R 1 9Q 
l \j ±ft±u — ±zy 


1 9Q 
u. -Lzy 


6944 


± 


415 


Fl 8 LI 


44.74 


± 


0.12 


8.40 


± 


0.12 


PP 1 49R4-01 ^ 


08fi 

U.UOO 


4976 


± 


731 


Fl 8 LI 


45.21 


± 


0.10 


8.44 


± 


0.15 


PP 1 4404-3^fi 


0.077 


2913 


± 


295 




44.75 


± 


0.07 


7.65 


± 


0.11 


PG 1444+406 


0.267 


4226 


± 


425 


12 


45.58 


± 


(0.10) 


8.56 


± 


0.12 


PG 1501+106 


0.036 


4673 


± 


470 


1 


44.22 


± 


0.09 


7.69 


± 


0.12 


PG 1512+370 


0.371 


8333 


± 


835 


1 


45.48 


± 


0.07 


9.08 


± 


0.11 


PG 1534+580 


0.030 


4987 


± 


500 


1 


43.86 


± 


0.10 


7.50 


± 


0.12 


PG 1545+210 


0.266 


4796 


± 


480 


F13,L14 


45.45 


± 


0.11 


8.58 


± 


0.13 


PG 1613+658 


0.129 


8073 


± 


30 


F1,8,L1 


45.09 


± 


0.07 


8.78 


± 


0.06 


PG 1704+608 


0.371 


3894 ± 


1469 


F1,8,L1 


45.56 


± 


0.18 


8.47 


± 


0.37 


PG 2130+099 


0.061 


3320 ± 


957 


F1,8,L1 


44.54 


± 


0.08 


7.62 


± 


0.26 


PG 2209+184 


0.070 


6595 ± 


660 


1 


44.41 


± 


0.17 


8.13 


± 


0.17 


PG 2251+113 


0.323 


3758 ± 


375 


F8,L14 


45.54 


± 


0.11 


8.43 


± 


0.13 


PG 2308+098 


0.432 


5328 ± 


535 


1 


45.87 


± 


0.08 


8.96 


± 


0.11 



a FWHM(C iv) measured in the single-epoch spectrum in units of km s 1 . 

b FWHM(Civ) and L A (1350A) are generally from the same reference, for the few exceptions: F: FWHM(Civ) reference, L: 
L A (1350A) reference. 

c log[AL A (1350A)/ergs s -1 ]. Errors in () are conservative errors assigned since no errors were quoted by the source study. 
d The central mass (and uncertainties; see text) estimated based on single-epoch UV spectroscopy, listed in logarithmic units. 

References. — (1) Wang, Lu, & Zhou 1998; (2) O'Brien, et al. 1998; (3) Koratkar, & Gaskell 1991; (4) B. Peterson, 2001, private 
communication; (5) Clavel, et al. 1991; (6) Korista, ct al. 1995; (7) Wanders, et al. 1997; (8) Wilkes, et al. 1999; (9) This author's 
own measurements of data by Vestergaard & Wilkes 2001; (10) Laor, et al. 1994; (11) This author's own measurements of published 
profiles by Marziani, et al. 1996; (12) Laor, ct al. 1995; (13) Wills, et al. 1995; (14) Schmidt, & Green 1983; Kellerman, et al. 1989. 



Black Hole Masses in Distant Active Galaxies. 



Table 3 

Line Width and Luminosity Regression Parameters 



Independent 


Dependent 


Sample 


N 


Slope 


Intercept 


Variable 


Variable 




(#) 


±error 


±error 


Bivariate Correlated Errors and Intrinsic Scatter Regression 


(BCES) 


Single-Epoch log L a 


Reverberation log L a 


A 


19 


1.03±0.05 


-1.59±2.60 


Bisector 




A 


19 


1.04±0.05 


-1.93±2.35 


Single-Epoch FW b 


Mean FW 


B 


18 


0.88±0.07 


172±181 


Bisector 




B 


18 


0.87±0.07 


206±187 


Single-Epoch FW 


Mean FW 


C 


17 


0.96±0.10 


~33±251 


Bisector 




C 


17 


0.94±0.09 


11±229 


Single-Epoch FW 


Rms FW 


c 


17 


0.97±0.09 


-237±266 


Bisector 




c 


17 


0.98±0.08 


-265±227 



a log [AL A (5100A)/ ergs s" 1 ] 

b FW = FWHM(H/3) ; Rms FW = FWHM(H/3, rms) ; Mean FW = FWHM(H/3, mean) 



Note. — Sample A: All PG quasars in Kaspi et al. (2000) 
Sample B: Sample A except PG1351+640 (no reverberation FWHM(H/3)) 
Sample C: Sample B except PG1613+658 (single-epoch FWHM(H/3) is not representative; 
see text) 



Table 4 

Mass Regression Parameters — Optical Measurements 

Independent Dependent Sam- N Slope Intercept 

Variable Variable pie (0) zberror ±error 

Bivariate Correlated Errors and Intrinsic Scatter Regression (BCES) 

logM BH (HS,S-E) a logM BH (H/3,rms) b (5 17 0.97±0.13 0.05±1.08 
BCES bisector C 17 1.03±0.14 -0.49±1.17 

logM BH (HS,S-E) logAf BH (H5,mean) b C 17 0.97±0.08 0.11±0.68 
BCES bisector C 17 0.91+0.11 0.59±0.91 

Sample C: All PG quasars of Kaspi et al. (2000) and MrkllO and Mrk335 (Wandcl 
et al. 1999) but excluding PG1351+640 [no reverberation FWHM(H/3) is available] 
and PG1613+658 [single-epoch FWHM(H /3) is not representative; see text] 



a M B H(H/3, S-E) = Single- Epoch mass estimates based on H/3 and optical continuum 
measurements. 

b log Mbh(H/3, rms) and log M B h (H/3, mean) are the reverberation masses based on the 
H P line widths only (in the rms and mean spectra, respectively) and the directly measured 
BLR sizes in Table 6 by Kaspi et al. (2000). 
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Table 5 

Mass Regression Parameters - UV Measurements 



Independent 
Variable 


Dependent 
Variable 


Sam- 
ple 


N 

(#) 


Slope 
±error 


Intercept 
±error 


Bivariate Correlated Errors and Intrinsic Scatter Regression (BCES) 


log[FWHM 2 (Civ) Ai A (K7 ] a 
BCES bisector 


log M BH b 


UVrev 
UVrev 


26 
26 


0.40±0.12 
0.79±0.17 


7.16±0.21 
6.52±0.36 


log[FWHM 2 (C IV) \L° x r ] 
BCES bisector 


logMeH 


UVrev,b 
UVrev,b 


25 
25 


0.54±0.19 
1.04±0.12 


6.91±0.33 
6.04±0.25 


log[FWHM 2 (Civ) XL\ h7 ] 
BCES bisector 


log M B h 


UV 
UV 


56 
56 


0.73±0.15 
1.03±0.16 


6.75±0.29 
6.20±0.32 


log[FWHM 2 (C IV) \L° x r ] 
BCES bisector 


logMeH 


UVb 
UVb 


55 
55 


0.91±0.12 
1.22±0.10 


6.39±0.22 
5.81±0.20 


logM BH ,uv(C IV) 


log Mbh (bisector) 


UVrev,b 


25 


1.07±0.14 


0.61±1.11 


logM BH ,uv(C IV) 


logMeH (bisector) 


UVb 


55 


1.23±0.10 


-1.78±0.81 



Sample UVrev: All Seyfert Is and PC quasars with reverberation masses from Wandel et al. (1999) 



and Kaspi et al. (2000) with published FWHM(C iv) and UV continuum luminosity 
Sample UVrev, b: Sample UVrev excluding NGC4151 

Sample UV: Sample UVrev supplemented with PC quasars from BG92 with published FWHM(C iv) 

and UV continuum luminosity 
Sample UVb: Sample UV excluding NGC4151 



a FWHM(Clv) is in units of 1000 kms 1 , and \L\is in units of 10 44 ergs s 1 and measured at 1350 A. 

b For the objects with available reverberation mass determinations log Mbh(H/3, rms) is used. For the 
additional PC quasars without such masses (sample UV and UVb only), the masses are estimated using the 
single-epoch Mbh(H/3); see section 5 



Table 6 

Probabilities of Mass Estimate Accuracies 



Calibration 


Factor 3 


Factor 6 


Factor 10 




(0.5 dex) 


(0.78 dex) 


(1.0 dex) 




accuracy 


accuracy 


accuracy 


Optical 


14/18 « 80% 


16/18 « 90% 


17/18 « 95% 


UV 


18/26 « 70% 


22/26 « 85% 


-23/26 « 90% 
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2000 4000 6000 8000 2000 4000 6000 8000 

Single-Epoch FWHM(HjS) [km s" 1 ] Single-Epoch FWHM(HjS) [km s" 1 ] 

Fig. 1. — The distribution of FWHM(H/3) in the mean (panel a) and rms (panel b) multi-epoch spectra and single-epoch FWHM(H/3) of 
the PG quasars presented by Kaspi et al. (2000) and MrkllO and Mrk335 (Wandel et al. 1999). The dotted line indicates a pure one-to-one 
relationship. The objects discussed in the text are labeled. The open squares denote Seyfert Is while triangles show measurements for the 
quasars. The three solid circles show the measurements of Boroson & Green (1992) which are based on spectra with both Fell emission and 
the narrow core component subtracted. These FWHM measurements are not always representative of the BLR velocity dispersion needed for 
this study (see text). The solid line is the best fit BCES bisector regression line based on all the objects in the diagram except PG1613+658 
(sample C; Table 3). The BCES (Y|X) and (X|Y) regressions are not plotted as they crowd the bisector. The single-epoch FWHM(H/3) 
scatter around a one-to-one relationship with FWHM(H j3, mean) to within 15% — 20% variation and around a similar relationship with 
FWHM(H j3, rms) to within 20% — 25%. Note, the ordinate range is different in the two diagrams. 
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Single Epoch log(AL A (5100A)/erg s" 1 ) Single Epoch log(AL x (5100A)/erg s" 1 ) 



Fig. 2. — (a) The distribution of the mean AL A (5100A) multi-epoch measurements (Wandel et al. 1999; Kaspi et al. 2000) based on 
monitoring data with respect to the single-epoch AL A (5100A) measurements of Neugebauer et al. (1987) and Schmidt & Green (1983). The 
errors in the mean AL A (5100A) are the rms around this mean. The errors in the sing le-epoch AL A (5100A) are propagated errors (see text). 
The short-dashed line (centrally positioned) denotes a one-to-one relationship. The dotted and long-dashed lines represent ±30% and ±60% 
luminosity variations, respectively (see text), (b) The best fit regression lines are shown for the BCES bisector (solid line). The BCES 
(Y|X) and (X|Y) regression lines would crowd the bisector, if plotted. All these BCES regression fits are consistent with a slope of 1.0. The 
single-epoch luminosities are offset by ±0.138 dex at AL A (5100A) Ri44.8 ergs s _1 , the mid-range luminosity, based on the BCES bisector. 




7 7.5 8 8.5 9 7 7.5 8 8.5 9 

log [M BH (Hj8, single-epoch) /M ] log [M BH (H/3, single-epoch) /M Q ] 

Fig. 3. — The reverberation masses derived from the rms (panel a) and mean (panel b) spectra plotted versus the single-epoch mass 
estimates based on optical spectral measurements. Triangles denote quasars, while squares denote Seyfert Is. The dashed line indicates a 
unity relationship. The dotted line is a BCES bisector regression line to all the PC quasars with H (3 measurements (shown), while the solid 
line is the BCES bisector when PG1613 is excluded (see text). As expected, Mbh(H/3, S-E) show less scatter with A/bh(H/3, mean). Both 
relationships are consistent with a one-to-one relationship to within the errors. 
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log [FWHM Z (CIV) (XL X (1350A)) - 7 ] -36.8 log [FWHM 2 (CIV) (AL X (1350A)) 07 ] -36.8 

Fig. 4. — Established and estimated central masses based on optical data plotted versus the UV measurements. This distribution is 
the basis of the calibration of the UV measurements. The "optical masses" (ordinate) in panel (a) consist of the reverberation masses, 
Mbh(H/3, rms), (filled, circled squares) derived from the rms spectrum and based on H/3 only. In panel (b) these masses are supplemented 
with the single-epoch mass estimates, Mbh(H/3, S-E), (filled squares) for the 30 PG quasars with no reverberation mapping masses. Dotted 
lines: BCES bisector regression lines to all the objects in each diagram. Solid lines: BCES(Y|X) and (X|Y) regressions to all the objects 
except NGC4151 (see text). Dashed lines: the BCES bisector for all objects excluding NGC4151. 




56789 56789 
log [M BH (CIV, calib) / M ] log [M BH (CIV, calib) / M ] 



Fig. 5. — The central mass estimates based on the calibration of the UV spectral measurements are compared to central masses measured 
(and/or estimated; panel b) based on optical data. Dashed line: pure unity relationship. Dotted lines: BCES(Y|X) and BCES(X|Y) regression 
lines. Solid line: the BCES bisector. NGC4151 was excluded from the regression analysis, (a) AfBH,uv estimates versus the established 
central masses, Mbh(H/9, rms) [UVrev sample only], (b) Mbh.uv estimates are plotted for the full UV sample versus the "optical masses" 
described in Figure 4b. For both diagrams the mass relationships are consistent with a unity relationship within the uncertainties. 
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Fig. 6. — The established central masses, Mbh(H/3, rms), based on optical multi-epoch spectral measurements plotted versus the deviations 
of mass estimates based on calibrated single-epoch spectra, (a) Deviations of masses estimated from optical spectra, i.e., A4"bh(H/3, S-E) 
divided by AfgH(H/3, rms). (b) Deviations in the masses based on UV spectral measurements, Mbh,uv(Civ). The uncertainties in the 
abscissa are the (propagated) uncertainties in the single-epoch masses (i.e., not the mass deviation error). A strictly unity relationship is 
indicated by the solid line. Offsets of ±0.5 dex (±ldex) are indicated by the dotted (dashed) lines. 
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Fig. 7. — The mass deviations from Figure 6 plotted here versus the H f3 and C IV line widths. The optical mass deviations are plotted 
against FWHM(H/3, rms), the width of the rms profile (i.e., the variable part), in (a) and versus the single-epoch FWHM(Civ) in (b). Lines 
and symbols are as in Figure 6. Notice that the larger mass discrepancies tend to occur for the most narrow lined objects in both cases (or 
those with the strongest variability occurring in the narrow line core: PG1613+658 and PG1704+608). 



